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Abstract 



This review summarizes recent developments in the understanding of high-energy cosmic 
rays. It focuses on galactic and presumably extragalactic particles in the energy range from 
^ <^ the knee (10^^ eV) up to the highest energies observed (> 10^^ eV). Emphasis is put on 

Qh observational results, their interpretation, and the global picture of cosmic rays that has 

Q emerged during the last decade. 
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1 Introduction 

Cosmic rays are ionized atomic nuclei reaching the Earth from outside the Solar System. Al- 
though already discovered in 1912, their sources and propagation mechanisms are still subject 
of intense research. During the last decade significant progress has been made and a consistent 
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picture of cosmic-ray observations begins to evolve. This review describes recent progress in the 
exploration of the sources and propagation of high-energy cosmic rays, focusing on observational 
results and the emerging global picture. 

Many reviews of cosmic-ray theory and observations are available in the literature. Most 
of these reviews concentrate on different aspects and energy ranges of high-energy cosmic rays 
which, in turn, are presented in much more detail than possible here. For example, covering the 
knee energy region, experimental data are compiled in [H [3] and a comparison of observations 
and model predictions can be found in [4J. The data in the range between the knee and the 
ankle, and their interpretation are discussed in [5j. Focusing on the upper end of the cosmic- 
ray spectrum, measurement techniques and observations are reviewed in El 1^- More 
theoretical and phenomenological aspects of the physics of ultra high-energy cosmic rays are 
subject of the reviews [H] and different source scenarios are discussed in depth in |18l [19j 
(acceleration scenarios) and [20] (non-acceleration scenarios). An exhaustive compilation of 
experimental results of the full cosmic-ray energy range can be found in \21j and a recent review, 
emphasizing measurement and analysis techniques, is given in [22j. 

In this article, we discuss high-energy cosmic-ray measurements covering the energy range 
from the knee to the highest energies. By concentrating mainly on observational results of the 
last decade and their implications for our overall understanding of cosmic rays, this review is 
complementary to the other articles. 

The exploration of cosmic rays is mainly driven by new experimental findings. Hence, we 
begin this review with a short historical overview, followed by an introduction to the physics of 
high-energy cosmic rays (Sec.[T]). 

In the energy region of interest cosmic rays are measured indirectly with large detector 
installations below the atmosphere, registering secondary particles produced in extensive air 
showers, initiated by high-energy cosmic rays. In Sec. [2} basic properties of air showers are 
introduced and major detection techniques are discussed. 

Recent experimental results concerning the flux of cosmic rays, their elemental composition, 
and studies of anisotropics in their arrival directions are presented in sections [3] to [5j The global 
picture evolving from these measurements and their impact on the present understanding of the 
origin of high-energy cosmic rays is emphasized in Sec.[6j The importance of the understanding 
of high-energy hadronic interactions for the interpretation of air shower data is underlined in 
Sec.[7| Concluding remarks (Sec.jsj) complete the review. 

Historical Overview 

Cosmic rays were discovered in the year 1912 by V.F. Hess during several ascends with hydrogen- 
filled balloons up to altitudes of 5 km [23j. He measured the ionization rate of air as function 
of altitude. Electrometers served as standard devices to measure ionizing radiation at this 
time . Hess found an increase of ionizing radiation with increasing height and he concluded 
that radiation penetrates from outer space into the atmosphere. For the discovery of the cosmic 
radiation, V.F. Hess has been awarded the Nobel Price in 1936. During the subsequent years W. 
Kolhorster made further ascends with improved electrometers, measuring the altitude variation 
of the ionization up to heights of 9 km [25j. 

In 1929 W. Bothe and W. Kolhorster measured coincident signals in two Geiger-Miiller 
counters [26j. Placing absorber material in between the two counters they also measured the 
absorption characteristics of the radiation. They concluded that the "Hohenstrahlung" (or 
cosmic radiation) is of corpuscular nature, i.e. consists of charged particles. Similar conclusions 
were drawn from measurements by J. Clay, who showed that the intensity of cosmic rays depends 
on the (magnetic) latitude of the observer [27j . This was a clear indication that a large fraction 
of cosmic radiation consists of charged particles. 
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Kolhorster continued his work with Geiger-Miiller tubes operated in coincidence. In February 
1938 he reported the discovery of coincident signals between two tubes with distances as far as 
75 m [28j. He concluded that the tubes were hit by secondary particles or showers generated 
by cosmic rays in the atmosphere. In the late 1930s P. Auger undertook investigations of 
cosmic radiation at the Jungfraujoch, Switzerland at 3500 m.a.s.l. He used Wilson chambers 
and Geiger-Miiller tubes separated by large distances and operated in coincidence [29]. Similar 
to Kolhorster, Auger concluded that the registered particles are secondaries generated in the 
atmosphere, originating from a single primary cosmic ray. 

In the 1940s the origin of the primary radiation could be revealed with measurements on 
balloons at high altitudes. M. Schein showed that the positively charged primary particles were 
mostly protons [30j. Cloud chambers and photographic plates were carried into the stratosphere 
and it was found that cosmic rays are made up of fully ionized atomic nuclei moving at speeds 
closely to that of light ^SJJ. Many nuclei of the periodic table up to Z ^ 40 were found and 
their relative abundances determined. Hydrogen and helium occur most frequently, and the 
distribution in mass of the heavier nuclei appeared to be similar to that in the solar system. 
Elements more massive than iron or nickel were found to be very rare. 

Since the mid 1940s large detector arrays were installed to measure extensive air showers. 
For most investigations detectors with a large surface and a short time resolution were required. 
Early detectors comprised Geiger-Miiller counters, progress in the development of photomultipli- 
ers lead to the application of scintillation counters and the newly available Cherenkov detectors. 
It was found that the energy spectrum of cosmic rays follows a power law dN/dE oc over 
a wide range in energy. In 1958 G.V. Kulikov and G.B. Khristiansen measured the integral 
electron number spectrum in air showers using an array of hodoscope counters [32j . They rec- 
ognized a kink in the spectrum around 6 x 10^ particles, corresponding to primary energies of 
several PeV (10^^ eV). This structure is now known as the "knee" in the energy spectrum. Since 
that time there is an ongoing debate about the origin of this structure. 

In the 1960s the air shower array of the M.I.T. group at Volcano Ranch, New Mexico was 
the largest cosmic-ray detector. The set-up comprised 20 stations equipped with scintillation 
counters, set up on a triangular grid, covering a total area of 12 km^. In 1962 the first event 
with an energy of about 10^^ eV has been recorded with this array [33j. Bigger air shower arrays 
were built (SUGAR [31], Haverah Park [35], Yakutsk [36], and AGASA [37]) and, after some 
initial attempts, the first successful fluorescence light detector, called Fly's Eye, was set up in 
Utah [38j. With these detectors, another feature of the cosmic-ray flux - first discussed in f39] 
- was firmly established in the early 1990s and is now known as the "ankle" H2]. 

Finally it should be mentioned that, in the early years, particle physics was done mainly by 
studying cosmic rays. In the 1930s investigations of the cosmic radiation lead to the discovery 
of new elementary particles such as the positron j43j or the muon [44j . The pion was discovered 
exposing nuclear emulsions to cosmic radiation at mountain altitudes in 1947 [45j. New unstable 
hadrons were found in cosmic-ray interactions in balloon-borne emulsion chambers [46] in 1971 
which were later, after the discovery of charm particles, identified as D mesons [U]- Also a 
number of exotic phenomena were observed [48j , none of which could be confirmed in accelerator 
experiments. 

Over time a standard description of cosmic rays evolved, which is briefly sketched in the 
following section. 

Nature and Origin of Cosmic Rays 

The solar system is permanently exposed to a flux of highly energetic ionized atomic nuclei - 
the cosmic rays. Their energies extend from the MeV range to at least 10^^ eV. The differential 
energy spectrum of all cosmic-ray particles is depicted in Fig.[T] It falls steeply as function of 
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Figure 1: All-particle energy spectrum of cosmic rays as measured directly with detec- 
tors above the atmosphere and with air shower detectors. At low energies, the flux of 
primary protons is shown. 



energy, decreasing by about a factor 500 per decade in energy. The flux decreases from more 
than 1000 particles per second and square meter at GeV energies to about one particle per m^ 
and year at a PeV, and further to less than one particle per km^ and century above 100 EeV. 

The strong decrease in flux poses a big experimental challenge and our knowledge about 
the particles and their origin is more and more limited with increasing energy (and decreasing 
flux). At sub-GeV energies individual isotopes are measured with small detectors in outer space 
and individual elements can be resolved with balloon-borne detectors in the TeV regime. At 
energies exceeding 100 TeV large detection areas are required to collect a suitable number of 
particles in a reasonable time. At present, such detectors are realized at ground level only 
and secondary particles generated in the atmosphere (the extensive air showers) are registered. 
At PeV energies, groups of elements could be resolved, while at the highest energies even a 
classification into "light" and "heavy" particles becomes already an experimental challenge. 

The energy spectrum follows a power law dN/dE oc over a wide energy range, indicating 
non-thermal acceleration processes. The spectrum is rather featureless as can be inferred from 
Fig.[l] However, small structures become clearly visible when the ordinate is multiplied with 
some power of the particle energy, as shall be discussed below (see e.g. Fig.[7|). The spectral 
index is 7 ^ —2.7 at energies up to several PeV. Then a steepening is observed, the so-called 
icnee, with 7 ^ —3.1 at higher energies. A further steepening, the second knee, occurs around 
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Nuclear charge number Z 

Figure 2: Abundance of elements in cosmic rays as function of their nuclear charge 
number Z at energies around 1 GeV/n, normalized to Si — 100 [50]. Abundance for 
nuclei with Z < 28 according to [5T]. Heavy nuclei as measured by ARIEL 6 [52j, [53j, 
HE AO 3 [51], SKYLAB [55], TIGER [56J, TREK/MIR |57], [58j, as well as UHCRE 
|59j . In addition, the abundance of elements in the solar system is shown according to 

m- 

4 X 10^^ eV. Finally, at about 4 x 10^^ eV, at the ankle^ the spectrum flattens again. 

The abundance of elements in cosmic rays is shown in Fig. [l] as function of the nuclear charge 
number. All elements of the periodic table have been found in cosmic rays. For the relatively 
abundant elements up to nickel, energy spectra for individual elements have been measured p4l3j. 
Abundances as obtained by several experiments at about 1 GeV/n are depicted. The cosmic- 
ray composition is compared to the abundance of elements in the solar system. Overall, both 
distributions look very alike. However, there exist certain differences, which reveal information 
on the acceleration and propagation of cosmic rays. 

The light elements lithium, beryllium, and boron as well as the elements below iron {Z = 26) 
and below lead {Z = 82) are more abundant in cosmic rays than in the solar system. They are 
assumed to be produced in spallation processes of the more abundant particles of the CNO, 
iron, and lead groups during the journey of cosmic rays through the Galaxy. Hence, they are 
frequently referred to as secondary cosmic rays. As the spallation cross section of the relevant 
nuclei is known at GeV energies, the ratio of secondary to primary cosmic rays is used to infer the 
propagation path length of cosmic rays in the Galaxy. An example is the boron-to-carbon ratio 
which has been measured as function of energy [61j. The ratio decreases as function of energy 
which is frequently explained in Leaky Box models by a rigidity-dependent [^decrease of the path 
length of cosmic rays in the Galaxy A(i?) = Aq{R/Ro)~^ . Typical values are Aq ^ 10 — 15 g/cm^, 

5 ^ 0.5 — 0.6, and Rq ^ 4 GV as reference rigidity. 

Cosmic-ray particles are assumed to propagate in a diffusive process through the Galaxy, 
being deflected many times by the randomly oriented magnetic fields (5^3 fiG). The nuclei are 
not confined to the galactic disc, they propagate in the galactic halo as well. The scale height of 
the halo has been estimated with measurements of the ^^Be/^Be-ratio by the ISOMAX detector 
[62j to be a few kpc. The abundance of radioactive nuclei in cosmic rays measured with the 
CRIS instrument yields a residence time in the Galaxy of about 15 x 10^ years for particles with 

^Rigidity is defined as particle momentum divided by its charge R[V] — z. 
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GeV energies [63j. 

The energy density of cosmic rays amounts to about per ^ 1 eV/cm^, a value comparable 
to the energy density of the visible star light psi ^ 0.3 eV/cm^, the galactic magnetic fields 
^ 0.25 eV/cm^, or the microwave background p^x ^ 0.25 eV/cm^. The power required 
to sustain a constant cosmic-ray intensity can be estimated as Lcr — PcrVl^esc ^ lO^-*^ erg/s, 
where Tg^c is the residence time of cosmic rays in a volume V (the Galaxy and the galactic halo). 
With a rate of about three supernovae per century in a typical Galaxy, the energy required 
could be provided by a small fraction 10%) of the kinetic energy released in supernovae. 
This had been realized already in 1934 by Baade and Zwicky [64J. The actual mechanism of 
acceleration remained mysterious until Fermi [65j proposed a process that involved interaction 
of particles with large-scale magnetic fields in the Galaxy. Eventually, this lead to the currently 
accepted model of cosmic-ray acceleration by the first-order Fermi mechanism that operates in 
strong shock fronts which are powered by supernova explosions and propagate from a supernova 
remnant (SNR) into the interstellar medium [66j. 

Diffusive, first-order shock acceleration works by virtue of the fact that particles gain an 
amount of energy AE" oc £" at each cycle, where a cycle consists of a particle passing from the 
upstream (unshocked) region to the downstream region and back. At each cycle, there is a 
probability that the particle is lost downstream and does not return to the shock. Higher energy 
particles are those that have remained longer in the vicinity of the shock and so have had time 
to achieve higher energy. After a time T the maximum energy achieved 

where Ps — Vs/c refers to the velocity of the shock. This results in an upper limit, assuming a 
minimal diffusion length equal to the Larmor radius of a particle of charge Ze in the magnetic 
fields B behind and ahead of the shock. Using typical values of Type II supernovae exploding 
in an average interstellar medium yields £'max ^ Z • 10^^ eV [[7^. More recent estimates give 
a maximum energy up to one order of magnitude larger for some types of supernovae \72\ . 
It has also been suggested that the cosmic rays itself interact with the magnetic fields in the 
acceleration region, yielding to an amplification of the fields, which in turn results in much 
higher energies that can be reached during the acceleration process [73]. With this mechanism 
cosmic rays are supposedly accelerated up to 10^^ eV. 

Information on the composition at the source can be obtained from measurements of the 
abundance of refractory nuclei. They appear to have undergone minimal elemental fractionation 
relative to one another. Comparing the derived abundance at the source to the abundance in 
the solar system reveals that the two samples exhibit a striking similarity over a wide range 
When uncertainties are taken into account, the abundances of particular isotopes are consistent 
with being within 20% of the solar values. This indicates that cosmic rays are accelerated out 
of a sample of well mixed interstellar matter. Hence, cosmic rays are "regular" matter, but 
accelerated to extremely high energies. 

The understanding of the origin of the knee in the energy spectrum is commonly thought to 
be a cornerstone for the understanding of the origin of (galactic) cosmic rays. Many approaches 
are discussed in the literature [1]. A popular explanation is that the knee is associated with the 
upper limit of acceleration by galactic supernovae, while the ankle is associated with the onset 
of an extragalactic population that is less intense but has a harder spectrum that dominates at 
sufficiently high energy. Another popular explanation is leakage of particles from the Galaxy. 
At energies in the GeV regime measurements indicate a decreasing path length of cosmic rays in 
the Galaxy. Extrapolating this to higher energies indicates that above a certain energy cosmic 
rays are not contained in the Galaxy (or the galactic halo) anymore. In a simple picture, this 
can be understood since the Larmor radius of a proton in the galactic magnetic field 
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Figure 3: Left panel: Energy loss lengths of protons and nuclei calculated for a redshift 
z = [86]. Right panel: Hillas plot [18j of astrophysical objects in which cosmic rays 
could be accelerated (adapted from [88]). 



becomes with increasing energy comparable to and finally exceeds the thickness of the galactic 
disk. 

If the knee is caused by the maximum energy attained during the acceleration process or 
due to leakage from the Galaxy the energy spectra for individual elements with charge Z would 
exhibit a cut-off (or a knee) at an energy — Z - Ec^ with the cut-off energy Ec for protons. 
The sum of the flux of all elements with their individual cut-offs makes up the all-particle 
spectrum. In this picture the knee in the all-particle spectrum is related to the cut-off for 
protons and the steeper spectrum above the knee is a consequence of the subsequent cut-offs for 
all elements, resulting in a relatively smooth spectrum above the knee. Since the abundance 
of ultra-heavy nuclei (see Fig.[T]) at GeV energies is very low as compared to iron the end of the 
galactic component is often assumed to be at energies around 30 x Ec. However, recently it has 
been proposed that ultra-heavy elements may play an important role at high energies [3J which 
yields a value of 92 x Ec for the end of the galactic component, coinciding with the second knee 
at 4 X 10^^ eV. 

Another interesting question is that of a natural end of the spectrum at high energy. Already 
40 years ago it has been realized that interactions of cosmic rays with photons of the cosmic 
microwave background (CMB) would result in a cut-off of the spectrum above 6x 10^^ eV All 
hadronic particles suffer significant energy losses during propagation above this energy. Protons 
interact with background photons forming mainly a A+(1232) resonance [79j and nuclei are 
broken up due to photodisintegration ^83j . A compilation of energy loss lengths of protons and 
various nuclei is shown in Fig. [3] (left) [8^. The energy loss length of photons depends on the flux 
of the universal radio background (URB) which is not well known (see, for example, discussion in 
[87j). Depending on the assumptions on the URB, the energy loss length is significantly smaller 
than or comparable to that of hadronic particles. 

The short energy loss lengths indicate that cosmic rays with energies above 10^^ eV should 
come from sources within a ~ 100 Mpc sphere. Astrophysical sources within our Galaxy are 
disfavored. Even though rapidly spinning young neutron stars could be thought of accelerating 
particles to the highest energies observed [89j, it would be difficult to explain the apparent 

^ Such a scenario has been pointed out first by Peters [75] and it has been suggested to call such a behavior a 
Peters cycle [76j. 
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isotropic arrival direction of UHECRs to energies beyond lO^*^ eV. As cosmic rays of energy 
greater than 10^^ eV are no longer confined by galactic magnetic fields, it is natural to assume 
that they are produced by extragalactic sources. 

Considering diffusive shock acceleration, which is thought to accelerate cosmic rays in SNRs, 
the magnetic field strength B in the source and the size R of the source region are related to 
the maximum acceleration energy by [18] 

£„„=.10«eVZA (2) 

where f3s is the shock velocity in units of c and Z is the particle charge. This relation is shown in 
Fig. [3] (right) for various astrophysical objects. The list of the very few viable candidate sources 
includes active galactic nuclei (AGN) [901 ED ES], radio lobes of FR II galaxies [93l [91], and 
gamma-ray bursts (GRBs) [95j (for a review of astrophysical sources, see [19j). 

Many alternative, non-acceleration scenarios - called top-down models - have been proposed. 
In these models, UHECRs are produced in decays of super-heavy objects such as super-heavy 
dark matter [98j, cryptons [lOOj, or topological defects All of these models postulate new 

particle physics and predict typically high gamma-ray fluxes at ultra-high energy |102| . Finally 
there are propagation model scenarios in which the GZK energy loss processes are evaded or 
shifted to higher energies. Examples are violation of Lorentz invariance [105J, the Z-burst model 
|110 ] or postulation of new particles with properties similar to protons |113^ I114| . Reviews of 
the different scenarios can be found in [20 1 flTSj . 

Measurements of the arrival direction distribution, primary mass composition and flux will 
be the key ingredients to solving the puzzle of UHECR sources. Given an expected angu- 
lar deflection of only a few degrees for particles above 10^^'^ eV in our Galaxy |116| and the 
existence of large cosmic voids with negligible magnetic fields [119J, high statistics measure- 
ments should finally allow cosmic-ray astronomy and reveal correlations with sources or source 
regions. Similarly, knowing the composition of UHECRs will restrict the classes of source mod- 
els. For example, a mixed composition would exclude top-down models. Another very impor- 
tant source of complementary information is given by secondary particle fluxes, i.e. gamma- 
rays and neutrinos, produced in UHECRs sources and during propagation (see, for example, 
[T22l[T23[[m ^ fT25l[T26l ). 



2 Detection Techniques 

Shower Properties 

In the following we will introduce some general properties of extensive air showers that are 
employed in cosmic-ray measurements. Detailed presentations of this subject can be found in 



Electromagnetic Showers 

Showers initiated by photons or electrons (called em. showers henceforth) are governed mainly 
by the particle production processes (i) bremsstrahlung of electron^ and (ii) pair production of 
electrons by photons. In addition to radiative energy losses, electrons are subject to ionization 
energy loss. The total energy loss dE/dX of electrons can be written as 

dE , . 



^In the following we will use the term electron to refer to both electrons and positrons. 
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where a{E) is the ionization energy loss given by the Bethe-Bloch formula which depends loga- 
rithmically on energy. The radiation length, Xr, depends on the material the shower evolves in 
and is Xji ^ 37g/cm^ in air. Particle multiplication and ionization energy loss are competing 
processes in showers. The energy at which an electron loses the same energy due to ionization 
and particle production is called critical energy, Ec = Xji{a{Ec)) ^ 86MeV. 

Some properties of em. showers can already be understood within the very simple model of 
Heitler [130] . In this model it is assumed that the incoming particle interacts in the atmosphere 
after traveling a depth distance Aem and produces two new particles with half the energy of the 
initial particle. These two new particles again interact at a distance Aem from their production 
point. After n generations of successive interactions the number of particles is 2^. The number 
of particles as function of depth X can be written as N{X) = 2^/^^™. The production of 
new particles continues until the energy of the secondary particles is smaller than the critical 
energy E^ below which absorption processes dominate over further particle multiplication. The 
maximum number of particles in the shower is A^max = Eq/Ec and the depth of shower maximum 
is given by X^ax = ^em ^0^2(^0 /Ec)^ with Eq being the primary particle energy. Heitler 's model 
is certainly an over-simplified picture of an air shower but illustrates two important features. The 
number of particles at shower maximum is approximately proportional to the primary energy 
and the depth of shower maximum increases logarithmically with energy. 

Approximate formulae for the longitudinal shower size profile and the lateral particle dis- 
tribution of em. showers have been derived from cascade theory |131t 1132] . Considering only 
shower particles of energy E, the depth of shower maximum is given by X^ax ^ ln(£'o/^)- 
Accounting for the energy distribution of particles in a shower this expression becomes (here 
written for electrons in a photon-induced shower) 



Amax ^ 



-'1)4 



(4) 



The increase of the depth of shower maximum per decade in energy is called elongation rate, 
DiQ. The elongation rate of em. showers is energy-independent, Df^^ = ln(10)Xi^ ^ 85g/cm^. 
For energies larger than E^ the energy spectrum of secondary particles in a shower follows 
approximately a power law dN/dE ~ £"-(1+^) with where s denotes the shower age parameter. 
The shower age is often phenomenologically defined as 5 3X/(A + 2 Amax)- Based on the 
solutions of the one-dimensional cascade equations given in |131| , Greisen developed a compact 
and still often used parametrization of the mean longitudinal shower size profile (number of 
charged particles) [ 134j 

A recent derivation of this expression is given in [13^J. In data analysis, a function proposed by 
Gaisser and Hillas |136j is often used to fit measured shower profiles 

The particle density in dependence on the distance to the shower core (lateral distance) is 
mainly determined by multiple Coulomb scattering of electrons. For electrons of low energy E, 
the increase of variance of the effective scattering angle per traversed depth 6X is approximately 
given by {69'^) ^ {Es/E)'^/XrSX, with Es ^ 21MeV [137J. Results of detailed calculations of 
the lateral particle distribution by Nishimura and Kamata [132J were parametrized by Greisen 
[T38j as 

=c(,mx)(^)'"(i + L)-'\ (7) 



rdrdip V^i/ \ 
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with C{s) = r(4.5 — 5)/[27rr^r(5)r(4.5 — 2^)] being a normalization constant. Eq. Q is called 
Nishimura-Kamata-Greisen (NKG) function. The lateral distribution of the shower particles 
scales with the Moliere unit ri = XrEs/Ec ^ 9.3g/cm^ and hence depends on the local air 
density. The effect of varying atmospheric density along the shower track can approximately be 
taken into account by calculating the Moliere unit not at observation hight but 2-3 radiation 
lengths higher up in the atmosphere |134j. 

At very high energy, two additional processes become important and change the characteris- 
tics of em. showers. First of all, the amplitudes of subsequent interactions of photons or electrons, 
which are independent at low energy, have to be added coherently [139j as the interaction length 
becomes comparable to the separation of subsequent interactions. The resulting effect is known 
as Landau-Pomeranchuk-Migdal (LPM) effect and leads to the suppression of new particle pro- 
duction in certain kinematic regions [141J. In air, the LPM effect typically becomes important 
at energies above 10^^ eV: Shower-to-shower fluctuations of em. showers increase drastically and 
the depth of maximum is shifted deeper into the atmosphere. 

Another important ultra high-energy process is that of magnetic pair production and brems- 
strahlung [143J. Gamma-rays of energies exceeding 10^^'^ eV can interact with the geomagnetic 
field of the Earth. Such interactions typically take place a thousand kilometers above the atmo- 
sphere. Mainly due to magnetic bremsstrahlung, a shower of more than 100 secondary photons 
and a few electrons is formed, which interact in the atmosphere simultaneously. Simulations 
of this effect can be found in \14:5\ I146^ 1147^ [148] . As the primary energy is shared by many 
secondary particles, the LPM effect hardly influences such showers. Due to the superposition of 
many lower-energy em. showers, shower-to-shower fluctuations of converted primary photons are 
significantly reduced. The dependence of the geomagnetic pre-shower effect on the local arrival 
direction can be used to search for photons in a model-independent way, see |145^ I148j . 



Hadron-initiated Showers 

Most of the differences between photon and hadron-initiated showers are related to the fact that 
hadronic showers develop a significant muon component whereas there are very few muons in 
purely em. showers. Furthermore, hadronic multiparticle production is characterized by large 
event-to-event fluctuations. 

Again a model similar to that from Heitler is useful to understand the basic properties of 
hadronic showers (see, for example, jl50( I15H 1152^ .153 J). For simplicity we assume that a 
hadronic interaction of a particle with energy E produces ntot new particles with energy E/ntoti 
two thirds of which being charged particles rich (charged pions) and one third being neutral 
particles rineut (neutral pions). Neutral particles decay immediately into em. particles (tt^ 27). 
Charged particles re-interact with air nuclei if their energy is greater than some typical decay 
energy E^^^ or decay otherwise. The number of generations of hadronic interactions, n, follows 
from £;dec = E^/{ntotT' 

Supposing that one muon is produced in the decay of each charged particle, one gets 

N,^M^^(^X , (8) 

V ^dec / 

with a — In rich/ In ntot ^ 0.86 . . . 0.93. The numerical values for a and E'dec depend on the muon 
energy threshold and are given in ||154| for different hadronic interaction models. The number 
of muons produced in an air shower increases almost linearly with primary energy and depends 
on the air density (through £'dec) ^ind the charged and total particle multiplicities of hadronic 
interactions. 

Also the energy transferred to the em. shower component can be estimated within this simple 
model. In each hadronic interaction, a certain fraction of the initial energy is transferred to the 
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em. shower component. After n generations the energy in the hadronic and em. components is 
given by 

Ehad ^ ( 3 ) = Eq - £^had- (9) 

Simulations show that the number of generations of charged pions is typicahy about 5 to 6 
|155j and increases shghtly with primary shower energy. Correspondingly the fraction of energy 
transferred to the em. component increases from about 70 — 80% at 10^^ eV to 90 — 95% at 
1020 eV. 

The depth of shower maximum of a hadron-induced shower is given by that of the em. shower 
component, X^^^. The first hadronic interaction produces em. particles of energy ^ E'o/ntot- 
Therefore one can write in lowest order approximation 



-^max (Eq) Ahad + -^max (^0 / ^t- 



otj 



where Ahad is the hadronic interaction length. From Eq. (10) follows the elongation rate theorem^ 
stating that hadronic showers have always an elongation rate equal to or smaller than that of 
em. showers [156J 

D'lo ^ {l-Bx-Bn)HlO)XRlniEo/Ec) 

^ {l-Bx-Bn)Dl^. (11) 

The coefficients are Bx = —dX^s^d/d^^E and Bn = dln{ntot)/dlnE (see also discussion in |154j 
for the next higher order term). 

If the primary particle is a nucleus, one can use the superposition model to deduce the main 
shower characteristics jl57j . In this model, a nucleus with mass A and energy £"0 is considered 
as A independent nucleons with energy = Eq/A. The superposition of the individual nucleon 
showers yields 

/Vmax (12) 

(13) 

^1-" • N^. (14) 

To a good approximation, while the number of charged particles at shower maximum is indepen- 
dent of the primary hadron, both, the number of muons and the depth of maximum do depend 
on the mass of the primary particle. The heavier the shower-initiating particle the more muons 
are expected for a given primary energy. In addition, the superposition of A independent show- 
ers naturally explains why the shower-to-shower fluctuations are smaller for shower initiated by 
nuclei as compared to proton showers. 

Detailed shower simulations conflrm the basic shower properties discussed so far. Quantita- 
tive predictions depend on the details of modeling particle production and transport [T58l[159]. 
Whereas em. interactions are rather well understood within perturbative QED, hadronic mul- 
tiparticle production cannot be calculated within QCD from flrst principles. Phenomenological 
models have to be used to describe the flnal states of hadronic interactions. The parameters 
of these models are determined by comparing the model predictions with accelerator measure- 
ments. Reflecting the different methods for describing data, low- and high-energy interaction 
models are distinguished. The former ones are typically based on the picture of intermediate 
resonance formation and decay as well as parametrizations of data. The latter ones are involving 
the production of color strings and their fragmentation. 
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Figure 4: Left panel: Longitudinal shower profile. Shown are the number of electrons 
with £^kin > 250 keV and muons with £^kin > 250 MeV (from jl51j ). Right panel: Mean 
depth of shower maximum calculated with different hadronic interaction models and 
CONEX (from [mj, modified). 



For high energy interactions (E'lab ^ 100 GeV), the hadronic interaction models often used in 
air shower simulations are DPMJET 11.55 and III ffGOlIIHI], EPOS [I62lll63], QGSJET 01 
and QGSJET II [1671, as weh as SIBYLL 2.1 [1571 [IMl [Eg . These models reproduce accelerator 
data reasonably well but are characterized by different extrapolations above a center-of-mass 
energy E'cms ~ 1.8 TeV {E ^ 10-^^ eV), leading to very different shower predictions at high 
energy [158. [T59l [1711 [T72] . 

The situation is different at low energy where more measurements from fixed target experi- 
ments are available. There, one of the main problems is the extrapolation of measurements to 
the very forward phase space region close to the beam direction and the lack of measurements of 
pion-induced interactions with light nuclei jl55j. At low energy, models based on data param- 
eterizations and/or microscopic models such as FLUKA [173J, GHEISHA [T7i] . UrQMD fT75] . 
or TARGET |176) are used. Differences of the predictions of these models are important for the 
number of muons [177j, in particular at large lateral distances. 

In Figs. [4][5} we illustrate some shower properties discussed only qualitatively so far by 
showing the results of detailed Monte Carlo simulations done with CORSIKA jl82| and CONEX 
|IS3j. 

The longitudinal profile of a typical proton shower of 10^^ eV is shown in Fig. [i] (left). In the 
region of the shower maximum, less than 1% of the charged shower particles are muons. The 
electromagnetic component is absorbed much faster than the muonic one. 

As muons are mainly produced in hadronic interactions, their number and lateral distribution 
can be used as composition-sensitive observables. The predicted muon distributions depend on 
the assumptions on hadron production in air showers. This is demonstrated in Fig. [5] where 
the expected number of muons and electrons is shown for showers initiated by proton, iron, 
and gamma-ray primaries, as calculated with different interaction models [172J. The energy 
dependence of the total muon number and ratio between the muon numbers of proton and iron 
showers is in agreement with the expectations from the superposition model. At high energy, 
the discrimination power of electron- muon number measurements is subject to large systematic 
uncertainties due to our limited understanding of hadronic multiparticle production. 

In Fig. [4] (right), we show the mean depth of shower maximum, (Xmax), for different primary 
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Figure 5: Expected number of muons and electrons in vertical showers at sea level. 
The curves show the FWHM of the distributions for different primary particles and 
energies, as obtained with QGSJET 01, QGSJET-II, and SIBYLL 2.1 [T72lfT85] . 



particles. The difference between proton and iron showers is in agreement with the predictions 
of the superposition model. Whereas the expected mean depth of maximum depends again sen- 
sitively on the chosen hadronic interaction model, the fluctuations are rather model-independent 
at a given energy. The elongation rate of em. showers obtained with detailed model simulations 
coincides with that predicted in em. cascade theory [131J. The increase of the em. elongation 
rate at high energy is caused by the LPM effect. Also, at ultra high-energy, photon interaction 
with the geomagnetic field can lead to a negative elongation rate. The elongation rates found 
for hadronic showers within different models (i?io^ ^ 50 . . .60g/cm^ at 10-*^^ eV) can be qual- 
itatively explained with the elongation rate theorem [156^ 1154] . An increase of the observed 
(-^max) faster than or comparable to the elongation rate of em. showers is a model-independent 
signature for a transition to a lighter composition. However, over a limited range in energy, 
drastic changes in the characteristics of hadronic multiparticle production can also lead to an 
elongation rate comparable to that of em. showers |186| . 

Finally, it should be mentioned that the em. shower component exhibits a number of uni- 
versality features that are independent of the primary hadron type and also rather insensitive 
to the primary energy. For example, near the shower core, the electron energy distribution is a 
universal function of shower age and the angular distribution of electrons depends almost exclu- 
sively on the electron energy and only slightly on shower age |189t 11901 1191| - Furthermore, the 
longitudinal shower profile (particle flux) at a given lateral distance can be well parametrized 
by an universal function of the depth of shower maximum and shower size at maximum jl92j . 
This reflects the universality of the lateral distribution of electromagnetic particles if the lateral 
distance is measured in Moliere units [194j. 

Measurement of Charged Secondary Particles 
Detection Techniques 

The classical set-up to measure air showers is an array of scintillation detectors, registering 
charged particles from the shower. In each detector the density of charged particles (mostly 
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electrons, positrons, but also some converted photons) is measured. From this information the 
lateral distribution of the electromagnetic component is inferred. This yields information on the 
position of the shower core and the total number of particles in the shower. Measurements of the 
arrival times of the particles with a resolution of a few ns allow to reconstruct the orientation 
of the shower plane and thus, perpendicular to it, the arrival direction of the primary cosmic 
ray. Due to the large number of secondary particles it is usually sufficient to cover only a small 
fraction of the total area with detectors. Typical values range from 1.2% for the KASCADE 
array to 5 x 10~^ for the Auger array. 

Examples for arrays in the knee region are the EAS-TOP experiment [ 196j . where 37 stations 
of scintillation detectors were distributed over an area of ^ 10^ m^, located above the Gran Sasso 
underground laboratory at an altitude of 2005 m, or the KASCADE experiment |197| . where 
252 detector stations are placed on a rectangular 13 m grid, covering an area of 200 x 200 m^. 
Detectors that operate at higher energy are the scintillator array of the KASCADE-Grande 
experiment [ 198| . or the AG AS A experiment, which covered an area of about 100 km^ with 111 
scintillation counter stations [37] . 

Another technique utilized to measure charged particles are water Cherenkov detectors, like 
the surface array of the Pierre Auger Observatory with its 1600 detectors distributed over an 
area of 3000 km^ |199| . They are relatively deep (typically ^ 1 m) compared to scintillators 
with typical thicknesses of the order of several cm. Consequently they have a larger response to 
showers with large zenith angles, e.g. the FWHM of the declination distributions were ^ 40° 
and ^ 75° for the Volcano Ranch (scintillators) [200J and Haverah Park (water Cherenkov) 
[35j experiments respectively. This results in a much larger sky coverage of a water Cherenkov 
experiment compared to a scintillator array. They are also suitable for the measurement of 
horizontal showers that can be used to detect neutrinos at the highest energies. 

Various techniques are applied for the detection of muons in air showers. Frequently, particle 
counters are covered with absorbers of lead, iron, or soil with a thickness of a few tens Xji to 
suppress the electromagnetic component, as in the KASCADE jl97| or AGASA [37] experiments. 
Muons can also be identified via their trajectories in tracking devices, like the HEGRA CRT 
detectors [201] or the muon tracking detector of KASCADE [202J. In underground laboratories, 
well shielded by rock, soil, water, or ice absorbers with a thickness corresponding to several 
1000 mw.e. (meter water equivalent), muons with thresholds in the TeV regime are registered 
[203] . 

The energy of hadrons is usually determined with calorimeters. The principle idea is to 
absorb an incoming particle and to measure the dissipated energy. Examples are the hadron 
calorimeters of the EAS-TOP [205] or KASCADE experiments [2(16] . 

Measured Parameters 

The direction of the shower axis and hence of the primary cosmic ray is obtained by measure- 
ments of the arrival time of the shower front in the detector plane. The shower direction is 
the normal to the reconstructed, slightly curved shower front. Typical detectors have a time 
resolution in the range from 0.5 ns up to a few ns. The angular resolution (68% value) of the 
KASCADE detector field decreases from 0.55° for small showers with Ne ^ 10^ to 0.1° for show- 
ers with Ne ^ 10^ [197J. The angular resolution of the Auger surface detectors decreases from 
values around 2.2° for showers with energies below 4 EeV to values from 0.5° to 1°, depending 
on the zenith angle, for higher energies (> 10 EeV) [207j . 

To obtain the position of the shower core and the number of particles, the particle densities 
are measured and an appropriate lateral distribution function is fitted to the data. The number 
of particles is calculated by integration of the measured lateral distribution. Historically, the 
choices of parameterizations of the electron and muon lateral distributions were influenced by a 
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review by Greisen [138j . The Nishimura-Kamata-Greisen (NKG) function became the standard 
function to describe the particle density p for a shower with the number of electrons and 
shower age s as function of distance r to the shower, see Eq. Q |132^ll38j . Traditionally, a fixed 
Moliere unit ri = 79 m is used and the parameter s is limited to 0.5 < s < 1.5. To parameterize 
the lateral distribution of muons, Greisen suggested the function [138J 

p,{r,N,) = CN,(^-Q + (15) 

to describe the density of muons as function of distance r to the shower axis for a shower 
with a total muon number A/"^. The Greisen radius is re = 320 m. 

The KASCADE group found that the lateral distributions of all three major shower compo- 
nents (electromagnetic, muonic, and hadronic) can be parameterized using the NKG function 
[208J. Fitting simultaneously the parameters N^, ri, and 5, the measured lateral distributions 
of the electromagnetic component can be reproduced with an accuracy of about 1%, yielding 
rf ^ 20 — 30 m. The experimental muon distributions are described with an accuracy of 5% 
using = 420 m. Finally, for the lateral distribution of hadrons with energies above 50 GeV a 
value Ti^lOm has been found. 

The NKG function has been analytically developed to describe pure electromagnetic showers. 
For hadron-induced air showers it exhibits shortcomings describing the measured electron lateral 
distributions, most obvious at large core distances. Investigations of the KASCADE group 
showed that the measured electron lateral distributions for showers with energies up to 10^^ eV 
and core distances up to 200 m can be described better using a modified NKG function [209J 

Pnkg = Ne ■ cis) .(^Y ^ (l + ^y ^ with c{s) ~ 



roj V ' roj 27rr^r(s - a + 2)r(a + /3 - 2s - 2) ■ 

(16) 

Optimizing the parameters with Monte Carlo data, the values a = 1.5 and /3 = 3.6 have been 
obtained, when ro = 40 m is used for the scale parameter. 

For showers in the EeV energy range, the AGASA group uses the parameterization 

to describe the lateral distribution of charged particles up to distances of several km from the 
shower axis j210j . inspired by a function suggested by Linsley [21 Ij. The parameters are 77 = 3.8, 
6 = 0.6 lb 0.1 and a Moliere unit tm = 91.6 m for near vertical showers with sec 9 < 1.2. 

Another parameterization for the electron lateral distribution has been suggested by the 
Haverah Park Collaboration |212j. It is also applied in the Auger Observatory [213J. The signal 
5 in a water Cherenkov detector is parameterized as function of distance r to the shower core 
as S{r) = A:r-(^+^/^-) for r < 800 m and S{r) = (l/800)^/cr-(^+^/^-)+^ at larger distances. 
The shape parameter 77 varies with zenith angle, while the parameter 6 and the scale radius 
Ts = 4000 m are fixed. 

The position of the shower core is determined by the KASCADE group with an uncertainty 
of 5 m for showers with an electron number ^ 10^ improving to less than 1 m for large 
showers {Ne ^ 10^-^) [197J. The error in the reconstructed number of electrons decreases in the 
same electron number interval from 18% to less than 4%. Similar values were obtained in the 
EAS-TOP experiment [196J. In the electron number range from = 10^*^ to 10^*^, the error 
in the position of the shower core improves from 7.5 to 2 m and the uncertainty in the number 
of electrons decreases from ANe/Ne = 28% to about 10%. 
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Energy Estimators 

One of the most important parameters to characterize a shower is the energy of the primary 
particle. Various methods are discussed in the hterature to obtain this value. 

For KASCADE, it has been found that, at sea level, the number of muons with energy above 
230 MeV in the range from 40 m to 200 m from the shower axis is a good measure for the 
primary energy independent of the mass of the primary particle [ 214^ . With its relatively high 
detector density the lateral distribution of muons is measured very well in the radius range from 
40 to 200 m. Extrapolating beyond these limits would introduce uncertainties related to the 
(less well known) shape of the lateral distribution. Hence, the number of muons is reconstructed 
using a distance range only, in which detectors are present. 

Another method is to use the correlation between the number of electrons and muons re- 
constructed. For example the CASA-MIA group uses the relation Eq = 0.8 GeV(A^e + 25A^^) to 
estimate the primary energy |215j . Similarly, the KASCADE-Grande experiment obtained the 
relation 



to estimate the primary energy as function of the observed number of electrons {Ee > 3 MeV) 
and muons {Ej^ > 300 MeV) at sea level for showers with zenith angle 9 j216j . 

A similar method has been applied by the AGASA group |217j . Here, the particle density 
as measured 600 m from the shower core is used to estimate the primary energy. At this 
distance the fluctuations in the lateral distribution are found to be relatively small [218^ I219| . 
The measured value is corrected for attenuation due to the zenith angle and converted to a 
value 50(600) for vertically incident showers. The conversion to primary energy via the relation 
E[eV] = 2.23 x 10^'^ • So(600)^-^'^[m^] yields at 10^° eV a difference between proton and iron 
induced showers of about 10% and a similar difference is obtained using two different codes to 
describe hadronic interactions in the atmosphere, namely QGSJET and SIBYLL [220] . 

The methods described so far depend on simulations of the shower development in the 
atmosphere. To avoid this uncertainty another method is being applied in the Auger Observatory 
j221 l|. It makes use of the constant intensity method, which relies on the fact that primary cosmic 
rays arrive isotropically. The value 5(1000) of the measured signal in a water Cherenkov detector 
at a distance of 1000 m from the shower axis is used. The dependence of this parameter on the 
depth in the atmosphere (which varies as sec 9 of the shower zenith angle 9) is obtained from the 
measured data. Using this dependence, the actually measured signal is converted to the value 
Sssi representing the signal for a shower with a zenith angle of 38°. Finally, the primary energy 
is estimated as £'[EeV] = 0.149 • S^^^ using the fluorescence telescopes as optical calorimeters, 
which define the energy scale in a direct way (see Sec. [2]) [221J. 

Composition Estimators 

To estimate the mass of the shower-inducing primary particle the following array observables 
are used: the electron-to-muon number ratio, the arrival time distribution of the particles, the 
curvature of the shower front, and the slope of the lateral distribution. 

The method applied most frequently is the measurement of the electron-to-muon ratio at 
ground level. Plotting the number of electrons and muons Njj^ in a plane as shown in Fig.[5j 
we find an energy axis (in the direction of the main diagonal) and (almost perpendicular to it) 
a mass axis. Using a Heitler cascade model the ratio of electrons to muons at shower maximum 
can be estimated, yielding the relation [153j 



logio 




0.313 logio + 0.666 log^g + 1.24/ cos 9 + 0.580 



(18) 




(19) 



This illustrates that logiQ Ne/N/j^ is a function of In A. It can be estimated that an uncertainty 
in logiQ{Ne /Njj,) of about 16% results in an uncertainty of about one unit in In. A. 

Also the shape of the shower front or the arrival time distribution of the particles at ground 
level is utilized as an indirect estimate of the depth of the shower maximum. For heavy nuclei, 
muons are produced earlier in the shower development and reach the ground also earlier as 
compared to the electromagnetic component whose particles follow the shower core and branch 
off to large lateral distances only very late before reaching the detector array. It is understood 
that the narrower pulse profiles correspond to higher production altitudes, which in turn lead 
to similar arrival times in the detectors (see e.g. [ 222^ i223j ). 

Measurement of Cherenkov Light 

Many particles in the shower disc travel with relativistic velocities through the atmosphere. 
Approximately one third of the charged particles emit Cherenkov light in the forward direction 
|224 j. the Cherenkov angle in air at sea level amounts to 1.3° only. Electrons (and positrons) 
are the most abundant charged species in air showers. Due to their relatively low Cherenkov 
threshold (21 MeV at sea level), they contribute mostly to the Cherenkov light in air showers. 

At present, for the detection of Cherenkov light, two techniques are applied: integrating 
detectors, in principle consisting of arrays of photomultipliers inside light collecting cones, look- 
ing upwards in the sky, and imaging detectors or telescopes, composed of large area collection 
mirrors and a camera with segmented read-out. Optical detectors such as Cherenkov detectors 
and fluorescence detectors (described in the next section) can only be operated during clear 
moon-less nights to obtain reliable data. This restricts their duty cycle to about 10%. 

Light Integrating Detectors 

The basic idea of integrating detectors is to measure the lateral distribution of the Cherenkov 
light with an array of photomultipliers distributed over a large area on ground level. To enlarge 
the collection area, the PMTs are installed inside light collecting cones (Winston cones). Such 
observations yield the lateral density distribution of Cherenkov photons. It can be parameterized 
by the empirical relation 



a combination of an exponential and a power law function [225]. C120 is the Cherenkov light 
intensity at 120 m distance from the shower core, s the exponential inner slope, and f3 the 
outer slope. The energy of the primary particle is strongly correlated with the photon density 
at 120 m, C120 grows approximately as E^-^^ . The average depth of the shower maximum 
-^max is approximately linearly related to the exponential slope s. Hence, from the Cherenkov 
measurements both, energy and mass of the primary particle can be derived. The latter through 
the dependence of the average depth of the shower maximum on the primary particle mass. 

This technique was pioneered by the AIROBICC experiment on La Palma island [226]. 
Another example is the BLANC A instrument, which was located at the Dugway Proving Ground 
in Utah, USA and operated in coincidence with the CAS A experiment at the same site |225| . 
BLANCA consisted of 144 angle-integrating detectors which recorded the lateral distribution of 
air shower Cherenkov light. A Cherenkov detector array is also installed in the Tunka valley 
close to lake Baikal in Siberia [2271 1229], consisting of 25 detectors that cover an area of 0.1 km2. 
It is planned to extend this installation with 133 optical detectors covering an area of about 




C120 • exp(s[120 - r/m]); 30 m< r < 120 m 
C120 • (r/120 m)-^; 120 m< r < 350 m. 



(20) 
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Figure 6: Left panel: Fluorescence light spectrum of air as measured by the AIRFLY 
experiment [2361 1237| with an electron beam of 3 MeV. The measurement was done 
at 20° C and a pressure of 800 hPa. Right panel: Ratio of calorimetric to total shower 
energy |184| . Shown are predictions calculated with CONEX |183| and different high- 
energy interaction models. The symbols present the mean values and la fluctuations 
are indicated by the shaded bands. 



Imaging Cherenkov Detectors 

Cosmic-ray events within the field of view of an imaging atmospheric Cherenkov telescope pro- 
duce a focal plane image which corresponds to the direction and intensity of Cherenkov light 
coming from the air shower. When the direction of the air shower core and the distance of the 
shower axis from the telescopes are known, simple geometry can be used to reconstruct the light 
received from each altitude of the shower. The amount of Cherenkov light produced is strongly 
correlated with the number of electrons in the shower and is used to estimate the shower size 
as a function of depth in the atmosphere from which the location of the shower maximum can 
be determined. This procedure is essentially geometrical and has the advantage of being almost 
independent of numerical simulations expect for the calculation of the angular distribution of 
Cherenkov light around the shower axis. 

Large Cherenkov telescopes are used to reconstruct air showers initiated by primary gamma 
rays in TeV 7-ray astronomy. The telescopes and analysis procedures are designed to effectively 
suppress the much more abundant (up to a factor of 1000) hadron induced showers. The 
technique was established by the pioneering work of the WHIPPLE telescope [231]. Among 
the presently largest instahations are the H.E.S.S. [232], MAGIC [233], and VERITAS [234] 
telescopes. 

An example for a Cherenkov telescope, optimized for the reconstruction of hadron induced 
showers was the Dual Imaging Cherenkov Experiment (DICE) [235J. It was located inside the 
CASA-MIA array in Dugway, Utah (USA) and comprised two telescopes, each equipped with a 
2 m diameter spherical mirror viewed by an array of 256 PMTs. 

Measurement of Fluorescence Light 

At very high energy {E > 10^^ eV), the fluorescence light technique can be used to measure 
directly the longitudinal profile of air showers. This technique is based on the detection of 
fluorescence light emitted by nitrogen molecules that are excited by charged particles traversing 
the atmosphere. 

There are two transitions of electronic states of the nitrogen molecule, called 2P and IN for 
historical reasons, that lead - in combination with the change of the vibrational and rotational 
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states of the molecule - to several fluorescence emission bands. A recent spectrometer measure- 
ment j236t I237| of these bands is shown in Fig. |6j The bands are labeled with the electronic 
transition type (2P or IN) and the change of the vibration quantum number [u^ v")^ for 
example, 2P(z/^, z/^'). Most of the fluorescence light emission is found in the wave length range 
between 300 and 400 nm (near UV). The lifetime of the excited states of nitrogen is of the order 
of 10 ns and the fluorescence light is emitted isotropically. 

The fluorescence yield, the number of photons produced per deposited energy, depends on 
the gas mixture in the atmosphere and atmospheric conditions. Collisions between molecules 
give rise to competing de-excitation processes (collisional quenching, see, for example, discussion 
in [238J). The importance of quenching increases with pressure and almost cancels the density 
dependence of the energy deposit per unit length of particle trajectory. This results in a weakly 
height dependent rate of about 4-5 fluorescence photons produced per meter and charged shower 
particle at relevant altitudes. During the last years several experiments have been set up for 
measuring the fluorescence yield under different atmospheric conditions j239^ I240| I24H 1242^ 
1236^243^1244^1245) . Recent progress is reviewed in j246| . Still, the fluorescence yield is currently 
known only to a precision of about 15% [240] as end-to-end calibrated data are not yet available. 

Recent measurements confirm the expectation that the fluorescence yield is independent of 
the energy of the exciting particle. Only at very low energy, deviations are expected [24T]. 
Therefore, if the atmospheric dependence of the fluorescence yield is taken into account, the 
fluorescence technique allows a calorimetric measurement of the energy deposited in the atmo- 
sphere. Simulations show that about 90% of the total shower energy is converted to ionization 
energy and, hence, is accessible for detection [248^ 1249^ [184 ]. The average ratio between the 
energy deposited in the atmosphere and the primary particle energy is shown in Fig. |6] (right). 
It depends on the primary particle type and energy and, to some extent, also on the model used 
in the simulations. However, as most of the shower energy is transferred to em. particles, the 
model dependence corresponds to an uncertainty of only a few percent of the total energy. In 
case of a gamma-ray as a primary, about 99% of the energy is deposited in the atmosphere. 

A complete reconstruction of a shower profile with the fluorescence technique requires the 
determination of the geometry of the shower axis, the determination of the Cherenkov light 
fraction, and the correction for the wavelength dependent atmospheric absorption of light. 

In shower observations with one fluorescence telescope (monocular observation), a shower- 
detector plane is given by the image of the light track. The orientation of the shower within this 
plane follows from the time sequence of the PMT signals [38 1 1250] . The angular uncertainty of 
the orientation of the shower-detector plane depends on the resolution of the fluorescence camera 
and the length of the measured track. Typically a resolution of the order of 1° is obtained. In 
general, the reconstruction resolution of the angle within the shower-detector plane is much 
worse and varies between 4.5° and 15° (for example, see [251| ). The reconstruction accuracy 
can be improved considerably by measuring showers simultaneously with two telescopes (stereo 
observation). Showers observed in stereo mode can be reconstructed with an angular resolution 
of about 0.6° [251| . A similar reconstruction quality is achieved in hybrid experiments that use 
surface detectors to determine the arrival time of the shower front at ground [252^ I207j . 

In fluorescence measurements, the Cherenkov light signal of air showers can be considered 
as a highly asymmetric background contribution but can also be exploited as an independent 
signal [253| . Knowing the longitudinal shower profile, the Cherenkov light contribution to the 
detected signal can be estimated using parametrized electron energy distributions and models 
for the angular distribution of the emitted Cherenkov hght [38] 1190^ ll9H I254| . 

In general, fluorescence detectors require continuous monitoring of atmospheric conditions, 
in particular the measurement of the wavelength dependent Mie scattering length and detection 
of clouds [255t 1199^ |256J . The temporal variations of the density profile of the atmosphere can 
lead to additional systematic uncertainties of shower reconstruction, in particular of the depth 
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of shower maximum |257j . 

The first fully functional air shower fluorescence detector was the Fly's Eye experiment in 
Utah that began taking data in 1982 and was operated for 10 years [38]. The Fly's Eye detector 
was a setup of two stations, Fly's Eye I (67 spherical mirrors of 1.5 m diameter, viewed by 880 
PMTs in total) and Fly's Eye II (8 mirrors viewed by 120 PMTs). The Fly's Eye I detector 
had a total field of view (FoV) of 360° in azimuth and 90° in zenith. With the smaller FoV of 
about 90° in azimuth and an elevation range from 2° to 38° degrees, the Fly's Eye II station 
was designed to measure showers in coincidence with Fly's Eye I. In October 1991 the shower 
of the highest energy measured so far, E = (3.2 zb 0.9) x 10^^ eV, was detected with Fly's Eye I 
[258] . 

The successor to the Fly's Eye experiment, the High Resolution Fly's Eye (HiRes), took 
data from 1997 (HiRes I) and 1999 (HiRes II) to 2006. The largest data set of HiRes is that 
of monocular observations with HiRes I, a telescope consisting of a ring of mirrors with a FoV 
from 3° to 16° in elevation and full azimuth. HiRes II is built up of two rings of mirrors covering 
elevation angles up to 30°. The optical resolution of the HiRes detectors is with 1° x 1° per 
PMT much higher than that of Fly's Eye. 

The Pierre Auger Observatory combines the observation of fluorescence light using imaging 
telescopes with the measurement of particles reaching ground level in a "hybrid approach" [199J. 
The southern Auger Observatory (near Malargue, Argentina) is the world's largest air shower 
detector and comprises 1600 polyethylene tanks set up in an area covering 3000 km^. Each 
water Cherenkov detector has 3.6 m diameter and is 1.55 m high, containing 12 m^ of high- 
purity water. A radio system is used to provide communication between each station and a 
central data acquisition system. Four telescope systems overlook the surface array. A single 
telescope system comprises six telescopes, overlooking separate volumes of air. Each telescope 
has a camera with 440 PMT pixels, whose field of view is approximately 1.5°. One camera 
overlooks a total field of view of 30° azimuth x 28.6° elevation. 

In the northern hemisphere, the Telescope Array (TA) is located in Millard County, Utah, 
USA |259| . It covers an area of 860 km^ and comprises 576 scintillator stations and three 
fluorescence detector sites on a triangle with about 35 km separation, each equipped with twelve 
fluorescence telescopes. 

Measurement of Radio Emission from Air Showers 

An independent measurement technique to observe air showers is provided by means of detec- 
tion of radio- frequency electromagnetic waves emitted from showers. Coherent radio emission 
generated by extensive air showers was theoretically predicted by Askaryan in 1961 |261| and 
experimentally discovered by Jelly et al. in 1965 at a frequency of 44 MHz |262j . Over a period 
of time this phenomenon has been considered as an interesting alternative to traditional methods 
of detection of high-energy cosmic rays with energy CTeater than 10^'^eV. In the 1960s and 1970s 
the experimental and theoretical efforts in this direction had only limited success |263| . Modern 
experiments, such as CODALEMA |264| and LOPES [265J, aim at studies of radio emission 
from air showers using modern, improved instruments. First break-throughs have been achieved 
[266$ I267j . At present, also activities are under way to install prototype antenna systems at 
the site of the Pierre Auger Observatory in Argentina to investigate the possibility for radio 
detection of air showers at the highest energies [268]. For a review of recent developments, see 

In addition to experimental difficulties there remain questions concerning the quantitative ra- 
dio emission theory. Several mechanisms of radio emission generation in air have been identified 
after the pioneering work of Askaryan where the coherent Cherenkov radiation of the charge- 
excess was put forward [26 Ij. This radiation is very strong for showers developing in dense media 
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Figure 7: All-particle cosmic-ray energy spectrum as obtained by direct measurements 
above the atmosphere by the ATIC [Ml IM], PROTON [Ml, and RUNJOB [283] 
as well as results from air shower experiments. Shown are Tibet AS7 results obtained 
with SIBYLL 2.1 [285J, KASCADE data (interpreted with two hadronic interaction 
models) [286], preliminary KASCADE-Grande results [287], and Akeno data [2881 H2] . 
The measurements at high energy are represented by HiRes-MIA [289^ 1290] , HiRes I 
and II [291J, and Auger [221j. 



[270] . In the case of air showers there is also an alternative radiation due to the acceleration of 
charged shower particles in the Earth's magnetic field. It is called geosynchrotron mechanism 
and has been recently investigated in detail [273J. The interrelation between these two essential 
mechanisms is not clear at present. Hence, also combined efforts are in progress, performing 
accurate radio emission calculations within the framework of a unified approach [276J. 



3 Energy Spectra 

The all-particle energy spectrum extending from 10^^ eV up to the highest energies is shown in 
Fig.[7| The flux as obtained from direct measurements above the atmosphere (represented in the 
figure through results from ATIC, PROTON, and RUNJOB) extends smoothly to high energies 
in the air shower detection regime. The all-particle spectrum can be approximated by a broken 
power law oc with a spectral index 7 = —2.7 below Ej^ ^ 4 x 10^^ eV. At this energy, the 
knee, the spectral index changes to 7 ^ —3.1. 

In the following we consider in more detail two energy regions: galactic cosmic rays up to 
energies of about 10^^ to 10^^ eV and the extragalactic component at higher energies. 



^The exact energy of the transition from galactic to extragalactic cosmic rays is presently not known, however, 
it is generally assumed to be in the energy range indicated, see also Sec.l6| 
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Figure 8: All-particle energy spectra in the knee region. Results from direct mea- 
surements by Grigorov et al. [282J, JACEE [M], RUNJOB [2811, and SOKOL [293] 
as well as from the air shower experiments Akeno 1 km^ [2^, BASJE-MAS [294], 
BLANCA [225J, CASA-MIA [2T5], DICE [235J, EAS-TOP ^95j, HEGRA [226J, KAS- 
CADE electrons and muons interpreted with two hadronic interaction models [286J, 
hadrons [296] , and a neural network analysis combining different shower components 
[297] . MSU [298], Mt. Norikura [299J, Tibet AS7 [300] and AS7-III [285], as weh as 
Tunka-25 \W\. 



Galactic Cosmic Rays 

Many groups published results on the all-particle energy spectrum from indirect measurements 
in the knee region 10^^ eV). The spectra obtained are compiled in Fig.jsj The ordinate has 
been multiplied by E^'^. The individual measurements agree within a factor of two in the flux 
values and a similar shape can be recognized for all experiments with a knee at energies of 
about 4 PeV. Also shown are results for the all-particle flux as obtained by direct observations 
above the atmosphere approaching energies up to 1 PeV. In the region of overlap, the results 
from direct and indirect measurements are in reasonable agreement. Typical values for the 
systematic uncertainties of the absolute energy scale of air shower experiments are about 15 to 
25%. Renormalizing the energy scales of the individual experiments to match the all-particle 
spectrum obtained by direct measurements in the energy region up to almost a PeV requires 
correction factors of the order of 10% [3J. A remarkable result, indicating that behind an absorber 
of 11 hadronic interaction lengths or 30 radiation lengths the energy of the primary particle is 
determined with an absolute error of the order of 10%. One should keep in mind that the 
experiments investigate different air shower components, are situated at different atmospheric 
depths, and use different interaction models to interpret the observed data. Nevertheless, the 
systematic differences are relatively small and the all-particle spectrum seems to be well known. 

Up to about a 10^^ eV direct measurements have been performed with instruments above the 
atmosphere. As examples, results for primary protons, helium, and iron nuclei are compiled in 
Fig. [9) Recently, also indirect measurements of elemental groups became possible, as discussed 
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Figure 9: Energy spectra for elemental groups, caption see next page. 
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Figure [9j Cosmic-ray energy spectra for four groups of elements, from top to bottom: 
protons, helium, CNO group, and iron group. 

Protons: Results from direct measurements above the atmosphere by AMS [301J, 
ATIC [302J, BESS CAPRICE ^04\, HEAT M, IMAX [307], JACEE 

[308J . MASS [310], RUNJOB [281], RICH-II [MS], [3l3j, SOKOL [293], 

[314] , and fluxes obtained from indirect measurements by KASCADE electrons and 
muons for two hadronic interaction models [286] and single hadrons [315] . EAS-TOP 
(electrons and muons) [316] and single hadrons [317] . GRAPES-3 interpreted with two 
hadronic interaction models [318J, HEGRA [319], Mt. Chacaltaya [320J, Mts. Fuji and 
Kanbala [32l], Tibet burst detector (HD) [322] and AS7 (HD) [323]. 
Helium: Results from direct measurements above the atmosphere by ATIC [302], 
BESS CAPRICE [304], HEAT [305j, ^06J, IMAX [307J, JACEE [308J, MASS 

[309] . [3T0] . RICH-II [m], RUNJOB [281], [313J, SOKOL [293], [324J, and fluxes 
obtained from indirect measurements by KASCADE electrons and muons for two 
hadronic interaction models [286], GRAPES-3 interpreted with two hadronic interac- 
tion models [318J, Mts. Fuji and Kanbala [321J, and Tibet burst detector (HD) [322J. 
CNO group: Results from direct measurements above the atmosphere by ATIC 
(C+0) [325J, CRN (C+0) [326], TRACER (O) [327], JACEE (CNO) [328], RUNJOB 
(CNO) [2Mi^ SOKOL (CNO) [293], and fluxes obtained from indirect measurements 
by KASCADE electrons and muons [286J, GRAPES-3 [SBJ, the latter two give results 
for two hadronic interaction models, and EAS-TOP [316] . 

Iron: Results from direct measurements above the atmosphere by ATIC [325], CRN 
[326J, HEAO-3 [329J, [330J, [331J, TRACER [327J (single element resolution) and [332] . 
[306J, JACEE [292J, RUNJOB ^84j, SOKOL [293J (iron group), as well as fluxes 
from indirect measurements (iron group) by EAS-TOP [316J, KASCADE electrons 
and muons [286J, GRAPES-3 ^8J, and H.E.S.S. direct Cherenkov light [333]. The 
latter three experiments give results according to interpretations of the measured air- 
shower data with two hadronic interaction models, namely QGSJET and SIBYLL. 
The gray solid lines indicate spectra according to the poly-gonato model ||3j. 



below and the results of the KASCADE and EAS-TOP experiments are shown in the figures 
as well. Also results from other air shower experiments are shown. HEGRA used an imaging 
Cherenkov telescope system to derive the primary proton flux [319] . Spectra for protons and 
helium nuclei are obtained from emulsion chambers exposed at Mts. Fuji and Kanbala [321J. The 
Tibet group performs measurements with a burst detector as well as with emulsion chambers and 
an air shower array j322( 1334] . GRAPES-3 uses the correlation between the registered number 
of electrons and muons to derive energy spectra for mass groups [318J. The H.E.S.S. Cherenkov 
telescope system derived for the flrst time an energy spectrum measuring direct Cherenkov light 
[333]. This light is emitted by the primary nuclei in the atmosphere before its flrst interaction, 
i.e. before the air shower begins [335j . Results for iron nuclei are shown. 

Over the wide energy range depicted, the flux as obtained by direct measurements is smoothly 
continued to higher energies with the results of air shower measurements. Despite of the ex- 
perimental uncertainties and systematic differences between different experiments and different 
interpretations of air shower data using various air shower models, a clear picture of the spectra 
for elemental groups is evolving. It is evident that the knee in the all-particle spectrum is caused 
by a depression of the flux of light elements. The measurements follow power laws with a cut-off 
at high energies. The spectra according to the poly-gonato model are indicated in the flgures 
as lines. It can be recognized that the measured values are compatible with cut-offs at energies 
proportional to the nuclear charge Ez = ^ • 4.5 PeV [3l 1336] . 
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Figure 10: Integrated aperture of different high-energy detectors corresponding to the 



data shown in Fig. [TTJ The AGAS A aperture refers to ah air shower data with 9 < 45° 
up to May 2003. The two HiRes detectors have different data taking periods: HiRes I 
from June 1997 to June 2005 and HiRes II from December 1999 to August 2004 [291J. 
The HiRes stereo exposure is that used for the recent anisotropy study [337J. The 
Auger exposure of 7000 km^ sr yr refers to data taking during the construction from 
Jan 2004 until August 2007, excluding events close to the array boundary [221J. The 
integrated aperture of the Yakutsk array includes data taken from September 1974 to 
June 2001 for TIOOO, and September 1979 till June 2001 for T500 [338J. The exposure 
shown for SUGAR is based on the re-analysis of the 5 highest energy events reported 
in j339j and corresponds to 11 years of operation. The Fly's Eye exposure in stereo 
mode is taken from [340J. The integrated aperture of the data set used in [341] for 
calculating the Haverah Park flux is 7.39 x 10^^ m^ ssr. 



Extragalactic Cosmic Rays 

There are several detector installations that were designed to measure cosmic rays at the highest 
energies. At the high-energy end, the total aperture and observation time determine the statistics 
of expected events. At low energy, the acceptance range is given by the employed detection 
technique and typically the distance between the detectors. 

A compilation of the integrated aperture (i.e. total exposure) reported by experiments with 



data above 10^^ eV is shown in Fig. 10 At ultra high-energy, about two times more events are 
expected in the Auger data set than HiRes has collected in monocular or stereo mode. The 
Yakutsk experiment is expected to have slightly more than a fourth of the statistics of HiRes I 
and about one half of that of AGASA. 

The measurements of the cosmic-ray fluxes obtained with these exposures are shown in 
Fig. [IT] (error bars indicate statistical uncertainties only) . In case several analyses of the same 
data set are available, only the most recent results are included in the plot. The shaded area, 
depicting the results of the analysis of the Haverah Park data j341| . accounts for some systematic 
uncertainties by assuming extreme elemental compositions, either fully iron or proton dominated. 
The highest energy point (Fly's Eye monocular observation) corresponds to the highest energy 
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Figure 11: Comparison of flux measurements scaled by E^. Only statistical errors 
are shown. Shown are the data of AGASA ^TTl [3l2], Akeno |288j, Auger [221j, 
Fly's Eye [SlOlSOl, Haverah Park [331], HiRes-MIA [289l|290l, HiRes Fly's Eye [291] . 
MSU [S3], SUGAR [339], and Yakutsk [333]. Yakutsk T500 (trigger 500) refers to the 
smaller sub-array of the experiment with 500m detector spacing and TIOOO (trigger 
1000) to the array with 1000m detector distance. The data of the MSU array are 
included to show the connection of the high-energy measurements to lower energy 
data covering the knee of the cosmic-ray spectrum. 



event [258]. For sake of clarity upper limits are not shown. 

It is common to present the data multiplied by to enhance deviations from a dN/dE ~ 
E~^ power law. In this representation, characteristic features such as the second knee at about 
10^'^-^ eV and the ankle at 10^^-^ eV are better visible. However, it should be emphasized that 
scaling the flux with energy to some power (e.g. E^J{E) or E'^'^J{E)) is misleading and does 
not reflect correctly the uncertainties of the measurements. In this presentation, the statistical 
uncertainties cannot be separated from the systematic energy calibration uncertainties. The 
importance of the systematic uncertainty of the experimental energy scale is demonstrated in 
Fig. 12 in which fluxes with the nominal energy scale of the experiments are compared with that 
after a model-motivated energy shift has been applied [347' '345] . After shifting the energy scales 
of the experiments, a very good overall agreement of the different measurements is obtained. 
In particular, the ankle in the cosmic-ray spectrum is clearly seen. The good agreement is a 
non-trivial observation as the position of the ankle and the overall flux change in a correlated 
way in the ^'^-representation. 

The data sets of the HiRes and Auger measurements provide evidence for a flux suppression 
at ultra high-energy. The statistical significance of the flux suppression is difficult to specify 
unambiguously. If one compares to a power-law flux a significance of more than 5 a is found in 
each of the data sets [221'^ 291J. Both spectra can be well described by models with uniformly 
distributed sources and a GZK suppression [34811349 ]. A similarly good agreement between GZK 
model predictions and the Yakutsk data was shown by several authors (for example, j 347i 1350] ). 

Only the AGASA data seem to disfavor a flux suppression. Assuming uniformly distributed 
sources of UHE protons and treating the normalization of the expected energy spectrum as free 
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Figure 12: Flux of UHECRs as measured with the four detectors that have the largest 
exposures, namely Yakutsk [344J AGASA [2171 [3l2], Auger [221], and HiRes [291J. 
Left panel: Cosmic-ray spectra as derived by the Collaborations using the calibration 
of the detectors. Right panel: Cosmic-ray spectra after re-scaling of the energy scale of 
the experiments to obtain a common position of the dip, from [345^, "346] . The nominal 
energy scales of the experiments have been multiplied by 1.2, 1.0, 0.75, 0.625 for Auger, 
HiRes, AGASA, and Yakutsk, respectively. 



parameter AGASA expects to observe 1.8 events above 10 eV with 11 actually detected. This 
corresponds to a 4.5cr deviation from the GZK cutoff spectrum |342| . Other assumptions on 
the shape of the GZK proton spectrum lead to the prediction of 2.4 expected events [351] 1350) . 
corresponding to a deviation of 3.9cr from the GZK cutoff hypothesis. Due to uncertainties in 
the absolute energy scale and partially also low statistics, the differences between the spectra of 
the different experiments are of limited statistical significance only |352t I353J. 

It should be kept in mind that a sudden and drastic change of composition from hadrons 
to photons could give observational results similar to a suppression of the flux j 354j observed 
with fluorescence detectors, but not with surface arrays. So far there are no indications for very 
deeply penetrating showers that would be expected in such a case. 

Given the importance of the absolute energy assignment to a reconstructed shower we will 
summarize the current systematic uncertainties below. 

The energy reconstruction of showers detected with the AGASA array is based on the scin- 
tillator signal 5(600) at 600 m from the shower core, where shower-to-shower fluctuations are 
the smallest and the relation between the signal and the primary energy is almost composi- 
tion independent [219^ I355j . The systematic error of energy assignment is analyzed in [2T7] 
in detail, see also discussion [356] • AGASA finds a total systematic uncertainty of the energy 
assignment of about 18%. The main sources of uncertainty are related to shower phenomenol- 
ogy and the simulation of the relation of 5(600) to the primary particle energy. In particular, 
the observed discrepancy between the surface detector signal at 1000 m from the shower core 
and the fluorescence-based calorimetric energy measurement reported by Auger [357J indicates 
that currently available shower simulations do not allow to obtain an absolute energy scale with 
an systematic uncertainty smaller than 20%. Therefore it is not surprising that discrepan- 
cies between the experimentally observed attenuation length for 5(600) and that expected from 
simulations hamper a re-analysis of the AGASA data [359J . 

The shower energy determination applied in HiRes is based on the track length integral 
^cai = <^eff / ^{X) dX^ where N{X) is a fit to the shower profile using the Gaisser-Hillas function 
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([g]) and aeff denotes the mean ionization energy deposit [360J . With HiRes being a fluorescence 
detector, the energy reconstruction is closely related to properties of the atmosphere which is 
serving as calorimeter. At the same time, atmospheric properties also determine the aperture 
of the detector. The HiRes flux measurements (HiRes I and HiRes H mono) are found to 
have similar systematic uncertainties ^3611 1362] I291J. The main contributions to the systematic 
uncertainty of the energy scale are the absolute calibration of the PMTs (10%), the limited 
knowledge of the air fluorescence yield (6%), and atmospheric conditions (9%). About 10% 
uncertainty results from the rescaling of the measured calorimetric energy to obtain the total 
shower energy [360], see also Sec.[2j Adding the individual contributions in quadrature, the 
overall systematic uncertainty of the energy reconstruction amounts to 17% [362J. 

In contrast to surface arrays, the aperture of fluorescence detectors has to be determined 
by simulations. Sources of uncertainty are here varying atmospheric conditions, simulation of 
shower profiles and detector trigger thresholds, and the primary cosmic ray composition. The 
uncertainty due to varying atmospheric conditions, mainly that of the vertical aerosol optical 
depth (VAOD), has been estimated to contribute to the aperture uncertainty 15% |36HI362] . In 
a recent study the other, simulation-related sources of uncertainty were found not to contribute 
significantly to the overall flux uncertainty of 30% [363J . 

The technique employed in the Auger measurement of the flux combines the advantages of 
surface detector arrays with that of fluorescence detectors [22 Ij. The surface array operates with 
almost 100% duty cycle and the aperture can be calculated in a rather straight-forward way for 
energies well above the trigger threshold. Fluorescence telescopes allow the direct measurement 
of the calorimetric shower energy, however, their duty cycle is only about 13%. Using a set of 
well-reconstructed hybrid eventsj^ one can calibrate the energy estimator for surface detector 
data in an almost model-independent way. This is done in a two-step process. First the shower 
signal at 1000 m, 5(1000), is corrected for attenuation to that of an equivalent shower of 38° 
zenith angle, Sss- To avoid any possible bias from simulations this is done with the constant 
intensity method by requiring the same number of showers per unit of sin^ 9. In the second step 
5^38 is converted to total shower energy. 

The aperture of the Auger detector increased continuously during construction and has 
an uncertainty of less than 3%. The systematic uncertainty of the energy scale coming from 
the fluorescence energy measurement is estimated to be 22%. The main contributions to this 
uncertainty are coming from the uncertainty of the fluorescence yield (14%), the calibration 
of the fluorescence telescopes (10%), and the reconstruction method (10%) [364" 221J. The 
described calibration procedure for relating S(IOOO) to the primary particle energy leads to an 
uncertainty of 7% at 10-*^^ eV increasing to 15% at 10^^ eV. 



4 Composition 

Mean Logarithmic Mass 

At energies below 10^"^ eV the abundance of individual elements has been measured with de- 
tectors above the atmosphere. At higher energies this is presently not possible due to the low 
flux values and the large fluctuations in the development of extensive air showers. Thus, in the 
past, mostly the mean mass has been investigated. An often-used quantity to characterize the 
composition is the mean logarithmic mass, defined as {In. A) = ^^r^lnA^, being the relative 
fraction of nuclei of mass Ai. Experimentally, (in A) is obtained applying two methods: (i) the 
quantity is proportional to the ratio of the number of electrons and muons registered at ground 



level {In A) oc logiQ{Ne/N^), see (19) and (ii) {In A) is proportional to the observed depth of 
the shower maximum, according to the relation X^^^ = -^max — XjilnA^ see (13). Hence, the 



^Events detected with both, the fluorescence telescopes and the surface detectors are called hybrid events. 
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Figure 13: Average depth of the shower maximum X^ax as function of primary energy 
as obtained by Auger [365], BLANC A |225], CACTI [366J, DICE [235], Fly's Eye 
[367 ]. Haverah Park [368], HEGRA [226], HiRes/MIA [290J, HiRes [369], Mt. Lian 
Wang [370J, SPASE/VULCAN [371J, Tunka-25 [229J, Yakutsk [372J. The hues indicate 
simulations for proton and iron induced showers using the CORSIKA code with the 
hadronic interaction model QGSJET 01 (— ), QGSJET H-3 (---), SIBYLL 2.1 (• • •), 
and EPOS 1.6 (— •)• 



maximum of an iron induced shower should be about 150 g/cm^ higher up in the atmosphere 
as compared to a proton induced shower (Xmax)- 



Recent measurements of the average depth of the shower maximum are compiled in Fig. 13 



covering energies from 10^ to almost 10^^ GeV. The experimental results are compared to pre- 
dictions of the average depth of the shower maximum from simulations for primary protons 
and iron nuclei. The CORSIKA code |182j has been used with the hadronic interaction model 
QGSJET 01 [161], QGSJET H-3 [167], SIBYLL 2.1 and EPOS 1.6 [373]. There are signif- 

icant differences between the predictions of the different models concerning the absolute values 
of Xmax- The differences become important when the model predictions are compared to the 
experimental data to derive information on the elemental composition of cosmic rays. 

Below 4x 10^ GeV the values obtained by different experiments exhibit a common trend, they 
seem to increase faster as function of energy than the simulations, which implies that the average 
composition would become lighter as function of energy. Above the knee (£'>4 x 10^ GeV) the 
measured values flatten up to about 4x10^ GeV, indicating an increase of the average mass in this 
energy range, as expected from sequential breaks in the energy spectra for individual elements, 
seen already in Fig.[9| Finally, above 4 x 10'' GeV the measured data exhibit about a constant 
slope for Xmax as function of energy. The slope is slightly steeper than the predicted slope for 
iron nuclei for all models shown. But, on the other hand, a comparison to the predicted proton 
values is not conclusive, while models like QGSJET 01 favor an extremely light composition at 
the highest energies, models like DPMJET 2.5 give a hint to an intermediate average mass. 

Knowing the average depth of the shower maximum for protons Xmax and iron nuclei X^^^ 
from simulations, the mean logarithmic mass is derived in the superposition model of air showers 
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Figure 14: Top panel: Mean logarithmic mass of cosmic rays derived from the av- 
erage depth of the shower maximum, see Fig.fTSl The hadronic interaction model 



QGSJET 01 is used to interpret the measurements. For comparison, results from 
direct measurements are shown as well from the JACEE |328j and RUNJOB |284| 
experiments. Bottom panel: Mean logarithmic mass of cosmic rays derived from the 
measurements of electrons, muons, and hadrons at ground level. Results are shown 
from CASA-MIA [371], Chacaltaya [375], EAS-TOP electrons and GeV muons [376]. 
EAS-TOP/MACRO (TeV muons) (STT], HEGRA CRT [20l], KASCADE electrons and 
muons interpreted with two hadronic interaction models [286] . hadrons and muons 
|378j . as well as an analysis combining different observables with a neural network 
[297] . and SPASE/AMANDA ^79\. 
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from the measured X^lf values using {In A) = - Xg,ax)/(X^L " ^max) • InAre- The 



corresponding {In A) values, obtained from the results shown in Fig. 13, are plotted in Fig. 14 



(top) as function primary energy using the hadronic interaction model QGSJET 01 to interpret 
the observed data. For comparison, also results of direct measurements are shown (JACEE and 
RUN JOB). 

In the figure three energy ranges may be distinguished for the indirect measurements. Below 
about 4 X 10^ GeV the individual experiments seem to indicate a decrease of {In A) with energy, 
while above this energy up to about 4 x 10^ GeV an increase with energy is exhibited. At the 
highest energies £^>4 x 10^ GeV, again a decrease with energy can be stated. 

Results of measurements of electrons, muons, and hadrons at ground level interpreted with 



the hadronic interaction code QGSJET 01 are compiled in Fig. 14 (bottom). They yield a 
clear increase of the mean logarithmic mass as function of energy. There seems to be some 
tension between the results obtained through the observation of the average depth of the shower 
maximum shown in the top panel and values derived from particle ratios measured at ground level 
depicted in the bottom panel [3l l380| . In particular, at energies below the knee (£" > 4 x 10^ GeV) 
the decrease of (in A) as derived from some X^ax measurements is not visible in the particle 
ratio results. 

Using the latest version of QGSJET II does not change the situation qualitatively. The three- 
fold structure of the results obtained is about the same as using QGSJET 01. The main difference 
are the absolute {In A) values which are shifted upwards by about 0.8 units for QGSJET II-3 
with respect to QGSJET 01. Using the hadronic interaction model SIBYLL 2.1 yields about the 
same values as for QGSJET II. Using lower inelastic hadronic cross sections in the QGSJET 01 
code and larger values for the elasticity of hadronic interactions the discrepancies between X^ax 
measurements and particle ratios at ground can be reduced |380^ 1381] . 



Spectra for Elemental Groups 

In addition to the mean mass as discussed above, it is interesting to investigate the energy 
spectra for individual elements or at least groups of elements. 

Information on the flux of primary protons can be inferred from the measurement of the 
flux of unaccompanied hadrons at ground level. With the K AS CADE hadron calorimeter the 
energy spectrum of single hadrons close to sea level has been measured in the energy range 
from 100 GeV up to 50 TeV [315j. Based on simulations using the CORSIKA code with the 
hadronic interaction model QGSJET 01 [164] the energy spectrum of primary hadrons in the 
energy range from 100 GeV to 1 PeV has been derived. Over the whole four decades in energy 
it can be described by a single power law. 

The KASCADE experiment used the number of electrons and muons (E"^ > 230 MeV) 
measured in the scintillator array to reconstruct energy spectra for five primary elemental groups 
[286J. Starting point of the analysis is the correlated frequency distribution of the number of 
electrons and the number of muons A/"^. Unfolding algorithms were applied to derive energy 
spectra for elemental groups. For the analysis the primary particles H, He, C, Si, and Fe were 
chosen as representatives for five mass groups. Details of the analysis and the used unfolding 
methods can be found in Ref. 1286 . 



The spectra obtained are presented in Fig. 15 To check the influence of the description of 



hadronic interactions in the atmosphere on the result, the same experimental data were unfolded 
using two interaction models, namely QGSJET and SIBYLL. The corresponding results are 
displayed in the figure. The resulting all-particle spectra for both models show a knee at about 
4 PeV and coincide within their statistical errors. The decrease of light elements across the knee, 
i.e. the occurrence of knee-like features in the light element spectra is revealed independently 



of the used simulation code, as can bee seen in Fig. 15 In contrast, the spectra of silicon 
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Figure 15: Cosmic-ray energy spectrum for five groups of elements as reconstructed by 
the KASCADE experiment using the hadronic interaction models QGSJET 01 (left) 
and SIBYLL 2.1 (right) to interpret the measured data [286J. 



and iron groups differ significantly and look quite unexpected. This can be understood by 
judging the ability of the simulations to describe the data. It turns out that both interaction 
models fail to reproduce the overall correlation between log^^g and log^^g observed in 

the data. In the case of QGSJET simulations the predictions are incompatible with the data 
in the low energy regime (simulations look too heavy), for SIBYLL incompatibility occurs at 
higher energies (simulations look too light). Summarizing this analysis, for the first time energy 
spectra for groups of elements were reconstructed from air shower data. The spectra indicate 
that the knee in the all-particle spectrum is due to fall-offs in the light element spectra resulting 
in a heavier composition above the knee. 

The EAS-TOP collaboration combined several detection systems to obtain information on 



the energy spectra of primary cosmic rays. The results are summarized in Fig. 16 



The proton spectrum in the energy range 0.5 - 500 TeV j317| has been obtained from the 
measurement of unaccompanied hadrons with a calorimeter, taking into account the contribution 
of helium nuclei as obtained by direct measurements |308|I382) . The proton spectrum is described 
over the whole energy range by a single power law. 

The helium and CNO fluxes in the energy region from 80 to 200 TeV have been studied from 
measurements of the Cherenkov light and TeV muons registered with the underground MACRO 
experiment [383J. Primaries are selected through their energy /nucleon by means of the TeV 
muon information. The shower energy is inferred from the measurement of the Cherenkov light 
yield at distances from 125 to 185 m from the shower core. The flux for p+He at 80 TeV 
and for p+He+CNO at 250 TeV was obtained. By subtracting the measured proton flux the 
following values were calculated: ^He{80 TeV) = (12.7 =b 4.4) x 10"'^ m~^sr~^s~^GeV~^ and 
$C7Vo(250 TeV) = (0.24 ±0.19) x 10"^ m-^sr-^s-^GeV-^. 

The all-particle energy spectrum is obtained from the measured shower size spectra in the 
knee region [295], showing the angular (i.e. depth) dependence of the knee position. The knee 
is observed at = 10^'^ in the vertical direction, corresponding to a primary energy Eq = 
(2 - 4) X 10^^ eV, and intensity of 10"'^ m ^s ^sr ^ with about 20% uncertainty. The obtained 
power law indices of the energy spectrum are: 71 = 2.76 =b 0.03 and 72 = 3.19 =b 0.06, below and 
above the knee, respectively. 

A composition analysis at knee energies was performed for vertical showers (1.00 < sec(O) < 
1.05) through measurements of the number of electrons and the number of muons with 
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Figure 16: Energy spectrum of cosmic rays as obtained by several detection methods 
applied in the EAS-TOP experiment [316J compared to results of direct measurements. 



energies above 1 GeV recorded in the muon detector at core distances r = 180 — 210 m (A^/iiso)- 
The experimental A/'^iso distributions, measured in different shower size intervals are fitted 
with simulated data to obtain energy spectra for groups of elements. Intrinsic fluctuations 
and measurement accuracies allow a three component analysis: light (constructed either with 
protons, and a mixture of 50% proton and 50% helium), intermediate (CNO), and heavy (Fe). 
The shaded areas in Fig. 16 indicate the energy dependence of the flux thus obtained for the 
three components. 

The Tibet air shower array has the advantage to be located at high altitude (4300 m, 
606 g/cm^). It comprises a scintillation counter array as well as emulsion chambers and burst 
detectors. The data were used to derive spectra for primary protons and helium nuclei, see 
Fig. [9] [322 \ 1334] . However, one has to keep in mind that only a few hundred events remain after 
quality cuts and are included in the analysis, which may indicate that the results are limited by 
their statistical significance. 

Similar to KASCADE and EAS-TOP, the GRAPES-3 experiment uses the correlation be- 
tween the registered number of electrons and muons to derive the flux for individual elemental 
groups. Spectra for protons, helium, as well as the CNO, silicon, and iron groups have been 
obtained [318J, see FigM 



Highest Energies 

In 1993 observations with Fly's Eye [40] gave the first indication of a systematic change of the 
cosmic-ray mass composition at very high energy. Analyzing the mean depth of shower maxi- 
mum, (Xmax), a change from an iron dominated composition at 10^^ eV to a proton dominated 
composition at 10^^'^ eV was found. However, an analysis entirely based on the mean X^ax is 
strongly model dependent (e.g. j3M IT58[ [5851 1386ll380l ). 

In Fig. [Tt] (left) a compilation of measurements of (Xmax) is shown together with model 
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Figure 17: Left panel: Compilation of fluorescence-based measurements of the mean 
-^max of very high-energy air showers. The data are from Auger [365], HiRes-MIA [289J, 
and HiRes stereo [369]. The model predictions are calculated with CORSIKA |182| 
and are taken from [386^ I185 j. The QGSJET predictions on the shower-to-shower 
fluctuations of the depth of maximum are indicated by the shaded (cross-hatched) 
area for proton (iron) primaries. Right panel: Fraction of iron, if data are interpreted 
with a hypothetical bi-modal composition of proton and iron primaries only. Shown 
are Akeno data [387], as well as results from HiRes [369J, Haverah Park [222J, and 
Vulcano Ranch [388^1389] . Upper limits to the iron fraction are obtained from AG AS A 
data [390] . 



predictions. Adopting the QGSJET 01 model |164j the conclusions of [IQ] still hold, though a 
mixed composition is expected at 10^^ eV. On the other hand, on the basis of models like SIBYLL 
2.1 [157^ 1169] 1170] or DPMJET 2.55 [160] a much more moderate change of the composition is 
derived. The model ambiguity of the interpretation of (Xmax) can be resolved to some degree 
by studying the measured distribution of X^ax [3841 1391j . 

In contrast to the old measurements of Fly's Eye [iQ] and Yakutsk |392| the HiRes data 
indicate a change from an iron-like to a proton dominated composition already at 10^^ eV. The 
two independent measurements are consistent in the overlap region. The large elongation rate 
of the low-energy data of ~ 93g/cm2 [210] can only be understood in terms of a change of 
composition (see Sec.|2]). Any model with scaling violations will predict a change to a lighter 
composition [156^ 1385^ I154 J. Also the muon density measured with the MIA detector |393j 
in the HiRes-MIA setup indicates a change from a heavy to light composition. The observed 
muon densities, however, are higher or similar than the expectation for iron primaries and not 



compatible with medium or light nuclei [290^ |2^ 

The new Auger data on the average depth of shower maximum [365] are, within the sys- 
tematic uncertainties of about 15g/cm^ in good agreement with the published HiRes data. On 
the other hand, if one would just analyze Auger data, a break in the elongation rate is found at 
a higher energy, at about 10^^*^^ eV. The elongation rate is 71 ± 5g/cm^ below this break and 
40 ± 4g/cm^ above [365]. The highest energy point of the Auger data is smaller than expected 
from a constant elongation rate. It is calculated from 13 observed showers and could indicate 
a transition to heavier elements as one would expect in case of a rigidity-dependent maximum 
energy of sources. 



The right panel of Fig. 17 shows a compilation of composition measurements based on particle 
detectors, as published in Ref. [390]. The Akeno result given here is re-scaled according to the 
predictions of the QGSJET model [359] Q 



^The original composition analysis of Akeno data [387] in the energy range 10^^"^ — 10^^"^ eV appeared to be 
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The muon density of ultra high-energy showers measured with AGASA was analyzed in 
Ref. g95]. Comparing the data with model predictions from AIRES ^396] and QGSJET 01 the 
following limits on a two component composition were derived: less than 35% iron in the energy 
range 10^^ - 10^^-^ eV and less than 76% at higher energy (90% c.l.) [39511390]. 

A re-analysis of Haverah Park data was done in Ref. [222 J. The authors use the sensitivity 
of the steepness of the lateral particle distribution to the shower development height, which 
in turn depends on the depth of shower maximum as composition sensitive observable. It was 
found that the predictions of CORSIKA with the QGSJET 98 model give a good description of 
the data if a two-component composition with about (66 ±2)% iron is used in the energy range 
from 2 X 10^'^ to 10^^ eV. At higher energy (from 10^^ to 2 X 10^^ eV), indications are seen for 
a transition to a lighter composition (see Fig. [l7|. This is supported by the number of inclined 
showers with E > 10^^ eV that have triggered the Haverah Park array j397] . The data analyzed 
in [397J agree well with simulations assuming all primary particles are protons, though no mass 
composition study was done. 

On the other hand, a first study of the time structure of Haverah Park showers with zenith 
angles less than 45° finds a more iron-dominated composition in the same energy range |223| . 
Also a re-analysis of Volcano Ranch data similar to the Haverah Park favors a large fraction of 
iron [398J. Using ^ 370 showers in the energy range from 5 x 10^^ to 10^^ eV a fraction of ~ 90% 
(75%) iron is found for a two component composition and QGSJET 98 (QGSJET 01). 

The discrepancy between muon density-based composition measurements and others based 
on features of the longitudinal profile underlines the shortcomings of the hadronic interaction 
models currently available. There seems to be a systematic deficit of muons predicted in simu- 
lations in comparison with data |290| 1357] . First progress in addressing this problem has been 
made in |163| by increasing the number of pair-produced baryons in the simulation. 



Ultra High-Energy Photons 

Using gamma-ray telescopes, such as MAGIC, H.E.S.S., or VERITAS, photons with energies 
of about 50 TeV have been observed from several sources in the sky, e.g. \4:W\ 1411] . At higher 
energies various air shower experiments searched for gamma rays. Air showers induced by 
primary photons develop an almost pure electromagnetic cascade. Experimentally they are 
identified by their relatively low muon content or their relatively deep shower maximum. Since 
mostly electromagnetic processes are involved in the shower development, the predictions are 
more reliable and don't suffer from uncertainties in hadronic interaction models. 

Results up to energies of 10^^ eV are summarized in Fig. 18 (left) [401J. There are no 



indications of a substantial fraction of gamma-rays in the high-energy cosmic-ray flux. 

The longitudinal profile of the highest energy event observed by Fly's Eye (£" ^ 3.2 x 
10^^ eV) [258] has been studied by several groups. Comparing the measured shower profile with 
Monte Carlo simulations shows that this event is well described by a hadronic primary particle 
|258| 1412] 1413] . However, due to the large reconstruction uncertainty of the atmospheric depth 
of the shower profile, a photon cannot be excluded j413| . 

The deeply penetrating muon component of inclined showers of hadronic origin is employed 
in an analysis of Haverah Park data in ^414^1397] . Using the primary cosmic-ray flux parametriza- 
tion of less than 48% of the observed events above 10-*^^ eV can be photons (95% c.l.). At 
energies above 4 x 10^^ eV this limit is 50%. 

Based on the analysis of muons observed in high-energy showers at AGASA the following 
upper limits were derived in [406^ »395] : 34%, 59% and 63% for primary energies above 10^^, 
lO-*^^-^^, and lO-*^^-^ eV, respectively (95% c.l.). A separate analysis of the 6 highest energy events 



in contradiction to the Fly's Eye composition interpretation "W, see also the analysis in [39] 
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Figure 18: Left panel: Upper limits for the photon flux derived from air shower ob- 
servations by the CASA-MIA [399J, EAS-Top [400J, and KASCADE [401j experiments 
compared to theoretical predictions |402j . Right panel: Upper limits on the fraction 
of photons in the integral cosmic-ray flux compared to predictions for GZK photons 
and top-down scenarios (SHDM, SHDM', TD) [403]. Experimental data are from the 
Auger surface detectors (Auger SD) [403] and a hybrid analysis (Auger FD) [404] , 
Haverah Park 005], AG AS A (marked AG AS A and AGASA') 006111071, AGASA and 
Yakutsk (A&Y) gOH], as weh as Yakutsk [409J. 



of AGASA was performed in |407| . Using a new method that accounts for the arrival direction 
of each individual shower, a limit of 67% at 95% c.l. could be derived for E > 1.25 x 10^^ eV. 

This new analysis method was also employed in a recent study of Auger shower longitudinal 
profile data [404J. No more than 16% photons are expected at 95% c.l. It should be noted 
that this limit as well as the one derived in |407j are independent of assumptions on hadronic 
multiparticle production at very high energy. It relies only on the simulation of photon showers 
which is much better under control than that of hadron-induced showers. 

In Ref. [408J, the scintillator signals of 10 showers above 10^^ eV from AGASA and Yakutsk 
were analyzed to derive a photon fraction limit. The authors report the strongest limit on the 
photon fraction for showers above 10^^ eV available so far, less than 36% photons at 95% c.l. 
To obtain this limit, new energies are assigned to the considered showers, partially being very 
much different from those originally reconstructed. This shows the importance of simulating 
and understanding detector effects which can only be done in a limited way in such studies. A 
confirmation by the AGASA and Yakutsk experiments would be very important to establish 
this limit beyond doubt. 

A compilation of recent upper limits on the contribution of photons to the all-particle flux is 



shown in Fig. 18 [403J. The best photon limits are the latest results of the Pierre Auger Obser- 
vatory [403] setting rather strong limits on the photon flux. They are based on measurements 
with the Auger surface detectors, taking into account observables sensitive to the longitudinal 
shower development, the signal rise time, and the curvature of the shower front. The photon 
fraction is smaller than 2%, 5.1%, and 31% above energies of 10^^ 2 • 10^^ and 4 • 10^^ eV, 
respectively with 95% confidence level. 

In top-down scenarios for high-energy cosmic rays, the particles are decay products of super- 
heavy objects. The decay process yields relatively high fluxes of photons, a typical feature of 
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Figure 19: Limits at 90% confidence level for a diffuse flux of assuming a 1:1:1 ratio 
of the three neutrino flavors at Earth [ 416^ [T3]. The experimental results are compared 
to predictions for GZK neutrinos and a top-down model [417J . 



such models |115| . Several predictions are shown in the figure [ 4T5t llOOj . These scenarios are 
strongly disfavored by the recent Auger results. It should also be noted that the upper limits 
are already relatively close to the fluxes expected for photons originating from the GZK effect 
^124j, shown in the figure as shaded area. 

Ultra High-Energy Neutrinos 

The detection of ultra high-energy cosmic neutrinos is a long-standing experimental challenge. 
Many experiments are searching for such neutrinos, and there are several ongoing efforts to 
construct dedicated experiments to detect them [418^ I4l9j . Their discovery would open a new 
window to the Universe [126J. However, so far no ultra high-energy neutrinos have been detected. 

m 

Due to interactions in the source region or during the propagation processes, ultra high- 
energy cosmic rays are expected to be accompanied by ultra high-energy neutrinos. The neu- 
trinos are produced with different abundances for the individual flavors, e.g. pion decay leads 
to a ratio \ — 2 \ \. However, due to neutrino oscillations the ratio expected at Earth is 
z/^ : z/^ : z/g = 1 : 1 : 1. 

To discriminate against the huge hadronic background in air shower detectors, neutrino can- 
didates are identified as nearly horizontal showers with a significant electromagnetic component 
|425j . The Pierre Auger Observatory is sensitive to Earth-skimming tau neutrinos that interact 
in the Earth's crust. Tau leptons from charged- current interactions can emerge and decay in 
the atmosphere to produce a nearly horizontal shower with a significant electromagnetic compo- 
nent. Recent results from the Pierre Auger Observatory together with upper limits from other 



experiments are presented in Fig. 19, Assuming an differential energy spectrum, the Auger 
Cohaboration derived a limit at 90% confidence level of EldNy^/dEy < 1.3 • 10"'' GeV cm~^ 
s"-*^ sr"-*^ in the energy range between 2 • 10^^ and 2 • 10^^ eV [416J. 

According to top-down models for ultra high-energy cosmic rays a large flux of ultra high- 
energy neutrinos is expected. As an example, the predictions of a model [41 7J are shown in the 
figure as well. This model is disfavored by the recent upper limits. It should also be noted that 
the current experiments are only about one order of magnitude away from predicted fluxes of 
GZK neutrinos. 



^Neutrinos produced in air showers (atmospheric neutrinos) [420 , in the Sun [4211 1422] , and during supernova 
1987A [4231 1424] have been detected, but are at energies weh below our focus. 
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Figure 20: Left panel: arrival direction of cosmic rays in celestial coordinates observed 
by the Super-Kamiokande experiment. Deviations from the average value for the same 
declination are shown (amplitude ±5%) |438| . Right panel: Cosmic-ray intensity map 
as observed by the Tibet experiment for various energy thresholds, from top to bottom: 
4 TeV, 6.2 TeV, 12 TeV, 50 TeV, 300 TeV pOj . 



5 Anisotropy 

The search for anisotropics in the arrival direction of cosmic rays on different angular scales 
can contribute to the understanding of the cosmic-ray origin, in particular the identification of 
source regions or individual sources. 

Galactic Cosmic Rays 

Large-scale anisotropics are connected to the propagation process of cosmic rays in the Galaxy, 
while small-scale anisotropics would be a hint towards cosmic-ray sources. However, it has to be 
considered that the Larmor radius ([T]) of protons with an energy around 1 PeV in the galactic 
magnetic field (5 = 3 /xG) is of the order of 0.4 pc. Hence, it is not expected to find any point 
sources for galactic cosmic rays. The situation changes for the highest energies, see below. 

An excess of charged particles with energies above 10^^ eV from the direction of a SNR 
(Monogem ring, d ^ 300 pc) has been reported |432| and later withdrawn |433| . This supernova 
remnant with an age of about 10^ yr has been suggested as possible single source of galactic 
cosmic rays [434J. A signal could not be confirmed neither by the KASCADE experiment [435] 
investigating cosmic rays with energies above 0.3 PeV nor by the Tibet experiment, looking for 
PeV 7-ray emission |436j . 

The Tibet air shower experiment has performed a northern sky survey, looking for TeV 7-ray 
point sources in a declination range from 0° to 60° [437J. A small excess is found, most likely 
caused by well known 7-ray sources such as in the Crab Nebula and Mrk 421. 

At higher energies, the KASCADE experiment has performed a detailed search for point 
sources, covering the whole visible sky at energies £"0 > 0.3 PeV j435j . The visible sky has been 
divided in cells with a size of 0.5°. Two distributions have been investigated, all events and a 
selection of muon poor showers. The latter have been investigated in order to look for potential 
gamma rays. They would manifest themselves in air showers with no or little muons only. Both 
investigations indicate an isotropic distribution of the arrival direction of cosmic rays. Special 
attention has been given to the region of the galactic plane, as well as to the vicinity of known 
SNRs and TeV- 7-ray sources. No significant excess could be found in either sample. 

The Super-Kamiokande experiment investigated large-scale anisotropics for cosmic rays with 
energies around 10 TeV j438j . The experiment registers muons from air showers with an energy 
threshold of 0.8 TeV. An anisotropy map in celestial coordinates is obtained by binning the 
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data in pixels of 10° x 10°. The result after applying a smoothing algorithm is shown in Fig. [20 



(left). Each cell represents the relative deviation from an isotropic distribution. A 3a excess 
("Taurus excess") is found with an amplitude of 1.04±0.20x 10~^ at right ascension a = 75°±7° 
and dechnation 6 = —5° =b 9°. On the other hand, a deficit ("Virgo deficit") is found with an 
amplitude of -(0.94 zb 0.14) x 10"^ at a = 205° ± 7° and 5 = 5° zb 10°. The angular difference 
between the Taurus excess and the Virgo deficit amounts to about 130°. A large-scale anisotropy 
would be expected if the Earth moves relative to a cosmic-ray rest system (Compton Getting 
effect [441]). For such a scenario the angular difference between maximum and minimum flux 
should be 180°. Within the relatively large angular uncertainties the anisotropy observed by the 
Super-Kamiokande experiment would be compatible with a Compton Getting effect caused by 
a velocity of about 50 km/s relative to the rest frame. This velocity is smaller than the rotation 
speed of the solar system around the galactic center (^ 200 km/s). This implies that the rest 
frame of cosmic rays (presumably the galactic magnetic fields) is co-rotating with the Galaxy. 

The Tibet experiment reported anisotropics in the same regions on the sky j440 j . The 
observed intensity for different energy thresholds is displayed in Fig. 20 (right) . For energies 
below 12 TeV the anisotropics show little dependence on energy, whereas above this energy 
anisotropics fade away, consistent with measurements of the KASCADE experiment in the 
energy range from 0.7 to 6 PeV |442| . A Compton Getting effect caused by the orbital motion 
of the solar system around the galactic center would cause an excess at a ^ 315°, 5 = 40° and 
a minimum at a = 135°, 6 = —49° with an amplitude of 0.35%. However, the measurements at 
300 TeV yield an anisotropy amplitude of 0.03% ±0.03%, consistent with an isotropic cosmic-ray 
intensity. Hence, a galactic Compton Getting effect can be excluded with a confidence level of 
about 5(7. This implies, similar to the result of the Super-Kamiokande experiment, that galactic 
cosmic rays co-rotate with the local galactic magnetic field environment. 

The Tibet experiment finds also an excess with a 0.1% increase of the cosmic ray intensity in 
the direction of the Cygnus region (a ^ 309°, 6 ^ 38°N) |440j . In the same region an excess of 
7-rays with energies around 10 TeV is seen by the Milagro experiment [443^ 1444] . Recently, the 
Milagro experiment reported also an excess of (charged) cosmic rays with a significance greater 
than 12a [445j . Presently, different explanations are discussed in the literature: the excess could 
be related to the Geminga pulsar as local cosmic-ray source [446J or cosmic rays from a local 
source could reach the Earth via a magnetic mirror [447]. 

At energies around the knee the Rayleigh formalism has been applied by several groups to 
characterize the large-scale anisotropy. The two-dimensional distribution of the arrival direc- 
tions of cosmic rays is reduced to one coordinate because of the limited field of view (typical 
experiments cover only a fraction of the whole sky) and the small amplitudes expected from the- 
ory. A first order approximation of the multipole expansion of the arrival directions of cosmic 
rays is a harmonic analysis of the right ascension values of extensive air showers. 

Applying the Rayleigh formalism to the right ascension distribution of extensive air showers 
measured by KASCADE yields no hints of anisotropy in the energy range from 0.7 to 6 PeV [442]. 
This accounts for all showers, as well as for subsets containing showers induced by predominantly 
light or heavy primary particles. Also other experiments, like Mt. Norikura [BE], EAS-TOP 
[339], Akeno [451J, and Adelaide [452J have derived Rayleigh amplitudes. Some experiments find 
anisotropics, however the phases do not agree between the different results. Hence, it seems to 
be more likely that all amplitudes derived should be considered as upper limits. 



Extragalactic Cosmic Rays 

Different groups have searched for large-scale anisotropics but no confirmed deviation from an 
isotropic arrival direction distribution has been found. A dipole amplitude that is compatible 
with full isotropy was found in an analysis of more than 135000 showers with energies from 
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3 X 10^^ to 10^^'^ eV of the Yakutsk array [453J . Also no significant anisotropy is found in the 
AGASA data |454| in the energy range between 10^^ to lO^^-^eV. 

In the energy range around 10^^ eV, however, AGASA found an excess of showers coming 
from directions near the galactic center and the Cygnus region [454] . Furthermore, analysis of 
SUGAR data [455J indicates an excess of cosmic rays coming from a similar region. The Haverah 
Park and Yakutsk arrays are located too far north to be able to see the excess regions of AGASA 
and SUGAR. 

With the galactic center being in the field of view of the southern Auger Observatory, it 
is ideally suited to search for these possible source regions, though the energy threshold for 
reaching 100% acceptance of the Auger surface array is just below 10^^"^ eV. Recently, the 
Auger Collaboration has performed an analysis of their low-energy showers from the direction 
of the galactic center. In the energy range from 10^^'^ eV to 10^^*^ eV, no abnormally over-dense 
regions were found in the neighborhood of the galactic center [456J. There are 506 (413.6) events 
found (expected) in the AGASA data set for a circle of 20° around 280° right ascension and —17° 
declination, corresponding to an excess ratio of 1.2. In the same region, the Auger Collaboration 
observed 2116 events while 2170 were expected. 

At energies above 10^^*^ eV, the large-scale structure of the arrival direction distribution 
appears, within the limited statistics of the AGASA array, isotropic [454J. This finding agrees 
with that of the HiRes Collaboration, performing a global anisotropy search based on ~ 1500 
events observed by HiRes I in monocular mode |457j . By combining data from arrays of the 
northern and southern hemispheres a full sky anisotropy study is done in [458J. Considering 
in total 99 showers from AGASA and SUGAR with E > 10^^*^ no large-scale anisotropy is 
found. The data set of the Auger Observatory analyzed in |459j contains more than 7000 events 
with E > 10^^-^ eV. No large scale anisotropy is found. The upper limit to the amplitude of a 
dipole-like anisotropy in right ascension is 3% with 95% CL ^459j. 

Small angle clustering could be an indication for point sources. At energies above 4 x 10^^ eV, 
clustering at an angular scale of 2.5° has been reported by the AGASA Collaboration [460J. This 
result could not be confirmed. An analysis of a HiRes I data set corresponding to an exposure 
similar to that of AGASA did not reveal any evidence for small scale clustering [463J. Also 
no significant clustering is seen in a data set from HiRes stereo observations of more than 
270 showers with E > 10^^ eV j464| and in the combined AGASA-HiRes stereo data set with 
£^ > 4 X 10^^ eV j465| . Clustering at the angular scale of the AGASA signal is also not found in 
the Auger data [466] . 

Uchihori et al. included showers above 4 x 10^^ eV from all four surface arrays of the northern 
hemisphere in their small scale correlation analysis j467j . The combined data set is found to 
contain many clusters, however, the statistical significance is low (^ 10%). On the other hand, 
there are indications for a correlation with the super galactic plane. Restricting the considered 
arrival directions to the range of ±10° off the super galactic plane the chance probability for 
finding doublets and triplets decreases to the order of 1%. 

Studies of small angle correlations with the Yakutsk array are difficult because of the angular 
resolution of the shower axis reconstruction, which is about 4° j468j . Nevertheless, at much lower 
energy, clusters of the arrival directions of showers in the energy range (1.3 — 4) x 10^^ eV were 
reported in [469j. The direction of these clusters seem to support a correlation with the super 
galactic plane. Dividing the observed cosmic-ray showers into isotropic and cluster components 
this correlation can be enhanced significantly [470J. 

A recent search for correlations at a medium angular scale has been carried out in [471J. 
Publicly available data sets from AGASA, HiRes, SUGAR, and Yakutsk with E > 4: x 10^^ eV 
were combined by adjusting their energy scales to shift the ankle in all data sets to the same 
energy. In addition events with E > 10^^ eV from Haverah Park, Volcano Ranch, and Fly's 
Eye were considered. In this set of 107 events in total, a signal at a scale of 25° is found in 
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Figure 21: Arrival directions (equatorial coordinates) of the highest energy cosmic rays 
observed with the Pierre Auger Observatory [484J (circles) and the HiRes telescopes 
j337j (squares). The asterisks indicate the position of active galactic nuclei (AGN) 
from the Veron-Cetty Veron catalog [481J up to a distance of 75 Mpc. The shaded 
area indicates the relative exposure of the Auger data set. The dotted line marks the 
galactic disk and the dashed curve is the super galactic plane. 

the autocorrelation function. Results from the Pierre Auger Observatory also indicate some 
clustering on intermediate scales with an angular separation of order of 15° to 25° [466J. The 
medium angle correlation reported in |466| has a 2% chance probability to originate from an 
isotropic distribution. 

There is a long history of searches for correlations with astrophysical point sources such as 
colliding galaxies and powerful radio galaxies. It appears almost impossible to assess unam- 
biguously the chance probability of such correlations. First of all, highly incomplete catalogs 
of astrophysical objects necessarily have to be used in these analyses. Secondly, the penalty 
factor for scans of several catalogs, selecting different classes of objects, distance and angular 
ranges, and other parameters cannot be calculated reliably. Part of these problems can be 
avoided by defining a prescription before analyzing a data set. This has been done by the Auger 
Cohaboration [172] . 

A correlation with BL Lacertae, at a distance exceeding the GZK energy loss length, has 
been claimed for the AGASA and Yakutsk high-energy data {E > 4: x 10^^ eV) |473j . The 
significance of this correlation is highly debated [476J as there might be "hidden" trials involved 
that cannot be corrected for with a Monte Carlo simulation. These correlations were tested 
with the independent data set of HiRes stereoscopic observations [480] . None of the previous 
claims could be confirmed. Recently, a ^ 0.8° correlation between BL Lacertae objects of the 
Veron-Cetty & Veron catalog ^481j with HiRes data {E > 10^^ eV) was pointed out [482J. This 
correlation was confirmed at a nominal significance of about 0.5% not only for the high-energy 
part but also for the entire set of HiRes stereo data |480) . An analysis of data from the Pierre 
Auger Observatory did not confirm a correlation of the arrival direction of cosmic rays with the 
positions of BL Lacertae objects in the southern hemisphere [483 J . 

A breakthrough in the anisotropy searches is the correlation discovered by the Auger Collab- 
oration (484|I485) . The arrival directions of the very highest energy events {E > 5.7 x 10-*^^ eV) are 
found to be correlated with the positions of active galactic nuclei (AGNs) from the Veron-Cetty 
& Veron catalog j481j . Out of 27 events observed with an integrated aperture of 9000 km^ sr yr. 
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20 are correlated with AGNs within an angular distance of 3.1°. The correlation was initially 
found in an exploratory search with different catalogs (12 out of 15 events were correlated). A 
prescription was set up to verify or reject the correlation hypothesis using an independent data 
set. Of the next 8 events that were detected, 6 were correlated with AGNs within the prescribed 
phase space, corresponding to a chance probability of incorrectly accepting the hypothesis of 
correlation of less than 1%. After accepting the correlation hypothesis, the Auger group made a 
scan to refine the correlation parameters and found an energy threshold of E^h = 5.7 x 10^^ eV, 
an "source" distance of less than z = 0.017 {D ^ 75Mpc) and a maximum angular difference of 
3.2° as optimal parameters. If one does not account for the penalty factor due to different trials 
and the parameter scans, the nominal chance probability of being compatible with isotropy of 
^ 10~^^ would be obtained. 

A sky map of the measured arrival direction distributions is shown in Fig. |2l] using equatorial 
coordinates. The relative exposure of the Auger Observatory is indicated by the shaded regions. 
The Veron catalog of AGNs is not a complete catalog. As expected, in the direction of the 
galactic plane the density of detected AGNs is lower than in other directions. 

The nearby AGNs are very good tracers of the local matter distribution. In particular, the 
super galactic plane is clearly visible. A recent update of the super galactic plane position even 
improves the correlation between AGNs and this matter over-density |486] . The correlation 
of the arrival direction distribution of UHECRs with nearby AGNs is also reflected by the 
autocorrelation function which shows some indications of anisotropy in the 15° to 25° range, as 
one would expect from the typical thickness of the super galactic plane. 

Given the limited statistics one cannot conclude from the found correlation that AGNs are 
sources of UHECRs. Subsequent studies of the published highest energy events of the Auger 
Observatory revealed correlations with the large scale structure of galaxies (487j. A correlation 
of UHECRs with the large scale structure in the cosmological neighborhood is also found in the 
Yakutsk data set [490J . 

The HiRes Collaboration has used the correlation parameters published by Auger to perform 
a search in their stereo data set. To obtain the same energy threshold as used in the Auger 
analysis, they scaled their reconstructed energies down by 10% j337j . Using all stereo data an 
exposure of roughly 4000 km^ sr yr is obtained at the highest energy (see Fig. 10). For a total of 
13 events two associations with AGNs were found, while 3.2 such correlations are expected for an 
isotropic arrival direction distribution. No correlation signal is found. Also the autocorrelation 
function of the highest energy events from HiRes is perfectly in agreement with the expectations 
for isotropy. The arrival directions of the selected 13 events are shown in Fig. 21_ as well. The 
exposure distribution of the HiRes data set is very similar to that of the Auger Observatory, 
but North exchanged with South. 

The discrepancy between the Auger and HiRes results are currently not understood, but it 
is clear that a possible difference of the energy scale of the two experiments could lead to such 
effects. In addition the reconstruction resolution has to be very good to reproduce the very 
sharp threshold of the correlation found in Auger data. An independent data set of similar size 
as the published one will allow to test the anisotropy signal. 



6 Astrophysical Interpretation 

Galactic Cosmic Rays and the Knee 

The measurements indicate that the knee in the all-particle energy spectrum is caused by a 
break in the spectra for the light elements, yielding an increase of the mean mass of cosmic rays 
in this energy region. Several scenarios are discussed in the literature as possible origin for the 
knee, see e.g. [4]. In the following, a current astrophysical picture of the origin of high-energy 
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Figure 22: Mean logarithmic mass as function of energy obtained by direct observations 
(dark gray area) and air shower experiments (hght gray area, obtained as weighted 
average of the results of many experiments, see ^) compared with different models 
(lines), a) Acceleration in SNRs (Berezhko et al. [491J, Kobayakawa et al. [492J, 
Stanev et al. 0931 SQl], Sveshnikova et al. j495| ): b) acceleration in GRBs (Plaga 
[i96] . Wick et al., gST] Dar gHE]), single source model (Eriykin & Wolfendale [499J), 
reacceleration in the galactic wind (Volk et al. |500j ): c) diffusion in the Galaxy 
(Kalmykov et al. |501j . Ogio et al. [502J . Roulet et al. j503| ): d) propagation in the 
Galaxy (Lagutin et al. |504j . Swordy [505j), as well as interaction with background 
photons (Tkaczyk [506J) and neutrinos (Dova et al. [507J). For details see [4J. 



cosmic rays is sketched, based on recent observations. 

One of the most popular explanations for the origin of the knee is that the spectra at the 
source exhibit a break. The bulk of cosmic rays is assumed to be accelerated in strong shock 
fronts of SNRs [66j. The finite lifetime of a shock front (^ 10^ a) limits the maximum energy 
attainable for particles with charge Z to £^max ^ ^ • (0.1 — 5) PeV. Many versions of this scenario 
have been discussed |49H 1493^ I494t I492( 1495] . The models differ in assumptions of properties of 
the SNRs like magnetic field strength, available energy, ambient medium, etc. The differences of 
the predicted (In A) can be inferred from Fig. 22 a. While older models [493J limit the maximum 
energy to about 0.1 PeV, recent ideas j495| . taking into account latest observations of SNRs, 
predict maximum energies above 1 PeV. In such a model sufficient energy is released from 
SNRs to explain the observed spectra. A special case of SNR acceleration is the single source 



43 



model [499j . which predicts in the knee region pronounced structures in the ah-particle energy 
spectrum, caused by a single SNR. Such structures can not be seen in the compilation of Fig.[8| 

In the literature also other acceleration mechanisms, like the acceleration of particles in 7-ray 
bursts, are discussed |496^I497^I498| . They differ in their interpretation of the origin for the knee. 
The approach by Plaga, assuming Fermi acceleration in a "cannon ball" is not compatible with 
the measured (in A) values, see Fig. 22 b. A different interpretation of acceleration in the cannon 
ball model yields - at the source - a cut-off for individual elements proportional to their mass 
due to effects of relativistic beaming in jets. The predictions of the actual model are compatible 
with recent data [498" "508]. However, it remains to be clarified how a detailed consideration 
of the propagation processes, e.g., in a diffusion model, effects the cut-off behavior observed 
at Earth. Gamma-ray bursts as a special case of supernova explosions are proposed [497J to 
accelerate cosmic rays from 0.1 PeV up to the highest energies (> 10^^ eV). In this approach 
the propagation of cosmic rays is taken into account and the knee is caused by leakage from the 
Galaxy leading to a rigidity dependent cut-off behavior. 

Also frequently discussed is the idea that the knee is a propagation effect. The propagation 
is accompanied by leakage of particles from the Galaxy. With increasing energy it becomes more 
and more difficult to confine the nuclei to the Galaxy. As mentioned above, the path length 
decreases as A oc E~^. Such a decrease will ultimately lead to a complete loss of the particles, 
with a rigidity dependent cut-off of the flux for individual elements. Many approaches have been 
undertaken to describe the propagation process [509^ I502| 1503^ 1505] I504j . The Leaky Box model 
[505] and the anomalous diffusion model [504] yield cut-offs significantly weaker than the data 
shown in Fig. |9] and [4J. 

The propagation as described in diffusion models jSOlj I502( I503| yields (ln74)-values which 



are presented in Fig. 22 c. The models are based on the same principal idea [509J, but take into 
account different assumptions on details of the propagation process, like the structure of galactic 
magnetic fields etc. This results in a more or less strong cut-off for the flux at the individual 
knees and, accordingly, in a more or less strong increase of (in A). The observed break of the 
spectra is relatively sharp. It seems to be difficult to generate such a behavior by a cut-off at 
the source or due to propagation only. More likely seems to be a combined approach assuming 
a break of the spectra at the source and leakage from the Galaxy, as e.g. discussed in 
During the propagation phase, reacceleration of particles has been suggested at shock fronts in 
the galactic wind |500| . Also this mechanism yields a rigidity dependent cut-off. 

Another hypothetical explanation for the knee are interactions of cosmic rays with back- 
ground particles like massive neutrinos [5071 l510j or photodisintegration in dense photon fields 
[506 , 1511] . Such models appear to be excluded with a high level of confidence. The interactions 
would produce a large amount of secondary protons, which results in a light mass composition 
at high energies, not observed by the experiments, see Fig. 22 d. Furthermore, a massive neu- 
trino, proposed in [507^ I510| can be excluded by measurements of the WMAP and 2dFGRS 
experiments j512j . 

A completely different reason for the knee is the idea to transfer energy in nucleon-nucleon 
interactions into particles, like gravitons |513j or extremely high-energy muons |514| , which are 
not observable (or not yet observed) in air shower experiments. The latter proposal seems to 
be excluded by measurements of the Baikal experiment [515J setting upper limits for the flux of 
muons above 10^ GeV. 

No point sources of charged cosmic rays were found in the knee region. The observed 
(large-scale) anisotropy amplitudes in the energy region of the knee are compared to model 
predictions in Fig. 23 Two versions of a leaky-box model [516], with and without reacceleration 
are shown. Leaky-box models, with their extremely steep decrease of the path length as function 
of energy (A oc E~^'^), yield relatively large anisotropics even at modest energies below 10^ GeV, 
which seem to be ruled out by the measurements, see also [ 516^ |4]. The measured values are 
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Figure 23: Rayleigh amplitudes as function of energy for various experiments, taken 
from 012]. The results obtained by Mt. Norikura gS], E AS-TOP [HH], Akeno [45 Ij, 
Adelaide [452] . and KASCADE |442| are compared to predictions of leaky-box models 
|516| and a diffusion model [517] . For the diffusion model, predictions for primary 
protons, iron nuclei, and all particles are displayed. 




almost an order of magnitude smaller. On the other hand, a diffusion model [517j, which is 
based on the idea of Ref. [509J, predicts relatively small values at low energies and a modest 
rise only. In the figure, predictions for pure protons and iron nuclei are shown together with 
calculations for a mixed composition. The predicted Rayleigh amplitudes are compatible with 
the measured values. This may indicate that diffusion models are a realistic description of 
cosmic-ray propagation in the Galaxy at PeV energies. 

In addition to information extracted from measurements of charged particles, important hints 
towards the origin of (hadronic) cosmic rays may be derived from observations of high-energy 7- 
rays. Recent observation of the H.E.S.S. experiment improve significantly the knowledge about 
galactic cosmic-ray sources. For the first time a spatially resolved image of a supernova remnant 
has been obtained with multi-TeV gamma rays as shown in Fig. [24] (left) |518j . The shell type 
supernova remnant RX J1713. 7-3946 has been studied in detail [520J. The image of the remnant 
has been divided into 14 regions to study the energy spectrum of gamma rays up to eight TeV. 
The indices of the power law spectra vary between 7 = 1.95=b0.08 and 2.24zb0.06 in the different 
fields with an overall index of the remnant 7 = 2.13 =b 0.03. The spectral indices are very close 
to the spectral steepness expected from Fermi acceleration at strong shocks. The morphology 
of the TeV 7-ray image agrees well with emissions in the 1-3 keV x-ray regime as measured by 
the ASCA satellite. 



The photon energy spectrum of the supernova remnant RX J1713 is presented in Fig. [24 
(right). Measurements in various energy ranges (ATCA at radio wavelengths, ASCA x-ray, 
EGRET GeV 7-ray, CANGAROO and H.E.S.S. TeV 7-ray) are compared to predictions of a 
model by Berezhko & Volk |519| . The solid line below 10^ eV indicates synchrotron emission 
from electrons ranging from radio frequencies to the x-ray regime. The observed synchrotron 
flux is used to adjust parameters in the model, which in turn, is used to predict the flux of TeV 
7-rays. The solid line above 10^ eV reflects the spectra of decaying neutral pions, generated 
in interactions of accelerated hadrons with material in the vicinity of the source (hadron + 
ISM ^ TT^ ^ 77). This process is clearly dominant over electromagnetic emission generated by 
inverse Compton effect and non-thermal bremsstrahlung, as can be inferred from the figure. The 
results are compatible with a nonlinear kinetic theory of cosmic-ray acceleration in supernova 
remnants and imply that this supernova remnant is an effective source of nuclear cosmic rays, 
where about 10% of the mechanical explosion energy are converted into nuclear cosmic rays 
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Figure 24: Left panel: 7-ray image of the supernova remnant RX J1713. 7-3946 ob- 
tained with the H.E.S.S. telescope. The superimposed contours show the x-ray surface 
brightness as seen by the ASCA satehite in the 1-3 keV range [518J. Right panel: 
Spatiahy integrated spectral energy distribution of the supernova remnant RX J 171 3. 
The solid line above 10^ eV corresponds to 7-ray emission from vr^-decay, whereas 
the dashed and dash-dotted curves indicate the inverse Compton and non-thermal 
bremsstrahlung emissions, respectively [519j . 



|519( 1521] . Further quantitative evidence for the acceleration of hadrons in supernova remnants 
is provided by measurements of the HEGRA experiment |522j of TeV 7-rays from the SNR 
Cassiopeia A |523| and by measurements of the H.E.S.S. experiment from the SNR "Vela Junior" 

m- 

In conclusion, it may be stated that a standard picture of the origin of galactic cosmic rays 
seems to emerge from the data. The measurements seem to be compatible with the assumption 
that (hadronic) cosmic rays are accelerated at strong shock fronts of supernova remnants. The 
particles propagate in a diffusive process through the Galaxy. As origin for the knee a combina- 
tion of the maximum energy attained in the acceleration process and leakage from the Galaxy 
seems to be favored. 



Transition Region 

Different scenarios are discussed in the literature for the transition from galactic to extragalactic 
cosmic rays. The transition most likely occurs at energies between 10^^ and 10^^ eV. 

The flux for elemental groups of the model of Hillas is shown in Fig. [25l [525] . The spectra are 
constructed with rigidity-dependent knee features at high energies. Reviewing the properties of 
cosmic rays accelerated in SNRs, and using the fluxes as derived by the KASCADE experiment 



(marked as component A in Fig. 25) Hillas finds that the obtained all-particle flux (dashed 
line, marked with Q) is not sufficient to explain the observed all-particle flux, see Fig. [25] [525]. 
Hillas proposes a second (galactic) component to explain the observed flux at energies above 
10^^ eV, marked as component "B" in the figure. An extragalactic component, marked as EGT, 
dominates the all-particle spectrum above 10^^ eV, for details see [525J. Very similar is the 
model proposed by Wibig and Wolfendale with a transition at higher energies between 10^^ and 
10-^^ eV [526) . In this model the galactic cosmic-ray flux extends to higher energies. Thus, a 
significant contribution of the extragalactic component is required beyond 10^^ eV only. 

Another possibility to match the measured all-particle flux is a significant contribution of 
ultra-heavy elements (heavier than iron) to the all-particle spectrum at energies around 4 x 
10^^ eV [3; 4J, as illustrated in Fig. 26 (left). The figure shows spectra for elemental groups with 



nuclear charge numbers as indicated, derived from direct and indirect measurements according 
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Figure 25: Breakdown of the cosmic-ray spectrum according to a model of Hillas [525] 
as the sum of galactic H, He, CNO, Ne-S, and Fe components with the same rigidity 
dependence, and extragalactic H + He having a spectrum oc E~'^'^ before suffering 
losses by cosmic microwave background and starlight interactions. The galactic com- 
ponents were given a turn-down shape based on a KASCADE knee shape as far as the 
point marked x. The dashed hne Q is the total galactic SNR flux if the extended tail 
(component B) of the galactic flux is omitted. [525j . 



to the poly-gonato model [3J. The sum of all elements is shown as solid line and is compared 
to the average experimental all-particle flux in the figure. In this approach the second knee is 
caused by the fall-off of the heaviest elements with Z up to 92. It is remarkable that the second 
knee occurs at E2nd ^ 92 • the latter being the energy of the first knee. In this scenario a 
significant extragalactic contribution is required at energies E^A x 10^^ eV. 

In the model of Berezinsky and collaborators [527^ I528j . the dip in the all-particle spectrum 
between 10^^ and 10^^ eV, see Fig. 26 (right), is interpreted as a structure caused by electron- 
positron pair production on cosmic microwave background photons p + 73^ ^ p + e+ + e~. 
Assuming a power law injection spectrum with a spectral index between 7 = —2.7 (without 
cosmological source evolution) and —2.4 (with cosmological source evolution), the spectrum can 
be described for E > 10^^'^ eV with a proton-dominated composition j527j . The shape of the 
dip is confirmed by data of the Akeno, AGASA, HiRes, Yakutsk, and Fly's Eye detectors after 
energy-rescaling |528| . Below a characteristic energy Ec ^ Ix 10^^ eV the spectrum flattens and 
the steeper galactic spectrum becomes dominant at < £'c. The transition energy Etr < Ec 
approximately coincides with the position of the second knee E2nd observed in the all-particle 
spectrum. The critical energy Ec is determined by the energy E^q = 2.3 x 10^^ eV, where 
adiabatic and pair-production energy losses are equal. Thus, the position of the second knee 
is explained in this scenario by proton energy losses on cosmic microwave background photons. 
The extragalactic component required in the poly-gonato model is somewhere between scenarios 
1 and 2 shown in Fig. [26] (right). It should be emphasized that the pair production mechanism 
requires the primary particles to be dominated (>80%) by protons [345 | 1125] . 

Traditionally, the ankle is interpreted as the characteristic signature for the transition be- 
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Energy Eg [GeV] e, Gev 

Figure 26: Left panel: Cosmic-ray energy spectra according to the poly-gonato model 
[3]. The spectra for groups of elements are labeled by their respective nuclear charge 
numbers. The sum of all elements yields the galactic all-particle spectrum ( — ) which 
is compared to the average measured flux. In addition, a hypothetical extragalactic 
component is shown to account for the observed all-particle flux ( — ). 
Right panel: Transition from galactic to extragalactic cosmic rays according to Berezin- 
sky et al. [527J. Calculated spectra of extragalactic protons (curves 1, 2, 3) and of 
galactic iron nuclei (curves 1', 2', 3') are compared with the all-particle spectrum from 
the Akeno and AGASA experiments. KASCADE data are shown as filled squares for 
the all-particle flux and as open circles for the flux of iron nuclei. 

tween galactic and extragalactic cosmic rays [525 1 1529] . In such a scenario, extragalactic cosmic 
rays dominate the flux above about 10^^ eV. This picture of the transition to extragalactic cosmic 
rays is supported by the pioneering observations of the Fly's Eye experiment that the compo- 
sition changes at about 10^^*^ eV [10 | I384j . New observations by HiRes-MIA and HiRes find a 
rather sharp transition from a heavy to a light composition at much lower energy, E ~ 10^^'^ eV. 
It is clear that the HiRes data are difficult to understand within a model in which naturally 
heavy elements should dominate the end of the spectrum of galactic cosmic rays just below 
10^^ eV. 

If one assumes that extragalactic cosmic rays are accelerated in processes qualitatively similar 
to those in our Galaxy then, at injection, the composition of extragalactic cosmic rays should 
be similar to that of cosmic rays at lower energy. Indeed, model calculations show that a mixed 
or even predominantly heavy source composition could, after taking propagation effects into 
account, be compatible with existing data [125( 1530] . 

On the other hand, the model by Berezinsky et al predicts a proton-dominated composition 
at energies as low as 10^^ eV. One of the advantages of this model is the natural explanation 
of the energy and the shape of the ankle. To obtain a good description of the ankle, there 
should not be more than ^ 20% He in the extragalactic cosmic-ray flux [125J 1345] . This could 
be interpreted as indication for either strong magnetic fields in the accelerating shock fronts or 
top-down source scenarios, which predict proton-dominated fluxes at not too high an energy. 

Understanding the nature of the ankle in the cosmic-ray spectrum has direct implications 
on the spectrum at much higher energy. For example, if the e^e~ pair production model is 
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confirmed one can conclude that (i) extragalactic cosmic rays are mainly protons, (ii) sources 
are cosmologically distributed, (iii) there should be a GZK suppression of the flux, (iv) an 
arrival direction correlation with local sources can be expected, and (v) the injection spectrum 
of extragalactic cosmic rays is rather steep {dN/dE ^ E~^;j > 2.4). 

Finally it should be noted that the neutrino flux contains complementary information for 
distinguishing between different model scenarios for the ankle [123^ [86^ l531| . 



Extragalactic Cosmic Rays 

Many authors assume extragalactic particles to be nuclei of intermediate to light mass. The 
discrepancy of elemental compositions derived from mean depth of shower maximum and elec- 
tron/muon number measurements, however, makes it impossible to use currently available com- 
position measurements as reliable criteria to disfavor models (see also discussion in ||533j). This 
is most strikingly seen in the prototype HiRes-MIA measurements [290^^ '289]. Whereas the 
mean depth of maximum data clearly showed the transition to a proton-dominated composition 
(QGSJET-based interpretation), the muon density at 600 m - of the same showers - appeared to 
correspond to primaries as heavy or even heavier than iron. A pure iron composition is obviously 
not compatible with data. The rather wide distribution of X^ax cannot be described with iron 
primaries only. 

The measurements of the HiRes experiment and the Pierre Auger Observatory have given 
evidence for a suppression of the flux at energies exceeding 4 • 10^^ eV p911|22l]. The question 
arises whether this steepening is due to the GZK effect, due to the maximum energy achieved 
during the acceleration processes, or due to both, a GZK energy loss process and an upper 
energy limit in the sources. 

The support for the existence of the GZK effect is provided by the correlation of the arrival 
directions with AGN, which imply a strong anisotropy of the arrival direction distribution. The 
anisotropy appears sharply at an energy of 57 EeV. At this energy, the flux measured by the 
Pierre Auger Observatory is about 50% lower than expected from a power law extrapolation from 
lower energies. Thus, there seems to be a connection between the steepening in the spectrum 
and the AGN correlation. However, if the observed suppression is due to the GZK effect one 
would expect either a light or rather heavy elemental composition above the GZK threshold. 
Intermediate mass nuclei are expected to break up very rapidly in interactions with the photons 
of the 3-K microwave background. The relative abundance of secondary protons produced during 
propagation according to a recent model [534| is displayed in Fig.[27| (right). Observations of 
the average depth of the shower maximum at the highest energies indicate a mixed composition, 
see e.g. Fig.[T3| 

The expected flux at Earth according to a propagation model of highest energy cosmic 
rays by Allard et al. |534| is presented in Fig.[27| (left). In this article different scenarios for 
the properties of the sources are discussed, like different elemental compositions and different 



spectral indices for the energy spectra at the sources. Figure [27| shows the result assuming a 
mixed elemental composition at the sources. The contributions of individual elemental groups to 
the all-particle spectrum are shown. The shape of the energy spectrum measured by the Pierre 
Auger Observatory is well represented by the all-particle spectrum derived from the model. 

Another important source of information are fluxes of particles produced during the propaga- 
tion of cosmic rays in the Universe. Fluxes of neutrinos and gamma-rays due to the interaction 
of ultra high-energy cosmic rays with the background photon radiation provide complementary 



information for discriminating models of UHECRs [ 122^ 1123^ I535|. I536| . In Fig. 19 measured 



^In the literature such particles are frequently called "secondary particles". They should not be confused with 
secondary particles generated in the Earth's atmosphere. 
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Figure 27: Results of a propagation model of highest energy cosmic rays by Allard et 
al. |534]. Left panel: Auger data compared to propagated spectra obtained, assuming 
a mixed source composition. The contributions of spectra for groups of elements to 
the all-particle spectrum are shown. Right panel: Expected relative abundance of sec- 
ondary protons as function of energy for different pure source composition hypotheses, 
assuming a spectral index of —2.1 at the sources. 

upper limits for the neutrino flux have been compared to a top-down model |417| . The latest 
Auger data disfavor this model. 

Photons of ultra high energy are easier to detect but their energy loss distance is very short. 
Nevertheless the measurement of the photon fraction in the primary cosmic-ray flux is one of 
the most promising methods of distinguishing between different source scenarios of extragalactic 
cosmic rays [2 mill5j . Recent experimental upper limits strongly disfavor predictions of top-down 
models j 415|ll00] . see Fig. [Ts] (right). Due to the down-cascading of photons in the extragalactic 
radiation background, fluxes of GeV photons are also a complementary source of information. 
For example, neutrino and photon fluxes of a Z-burst model are discussed in Ref. [122]. Both, the 
neutrino flux limit by the FORTE satellite [537J and the EGRET diffuse extragalactic photon 
flux analysis |538| severely constrain this model. Recently, also even more stringent limits were 
set by the ANITA experiment [539] , ruling out most of the parameter space of the Z-burst 
model. In summary, top-down models are disfavored by recent experimental results of various 
experiments. 

One of the key questions in the field of high-energy astroparticle physics is the understanding 
of the observed anisotropy above the GZK energy threshold. 

The angular scale of about 3° would favour protons if the AGNs within this correlation angle 
are indeed the sources of UHECRs [485J. This would be at variance with the Auger data on 
the mean depths of shower maximum. Most likely, the AGN correlation has to be interpreted 
merely as a signature of anisotropy and correlation with the nearby matter distribution |540| . 
First of all, the rate of misidentification of a potential source along the line of sight of the arrival 
direction of a cosmic ray is very large. Secondly, the AGNs within the correlation window are 
found to be often less powerful Seyfert 2 galaxies |541j. Thirdly, the AGNs in our cosmological 
neighborhood are strongly clustered, making it difficult to distinguish between a true AGN 
correlation and a correlation with larger scale structures. 

Instead of assuming a single source in each of the directions of the measured UHECRs, one 
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can assign most of the arrival directions to about three sources or source regions [542J. This 
model scenario would require rather weak galactic magnetic fields and particles in the mass 
range up to carbon. 

In all anisotropy studies even at the highest energies, the galactic magnetic fields play a 
central role [543J. Knowing their structure to sufficient detail would, for example, allow the 
determination of the charge of UHECRs. Other interesting applications are the search for the 
cosmological Compton-Getting effect, a 0.6% dipole anisotropy that is expected for cosmological 
sources [544J. 

A much higher number of UHECRs has to be collected for detailed anisotropy studies of 
the required sensitivity. If UHECRs are confirmed to be protons, the data can be used for 
proton astronomy including studies of energy spectra of individual sources and magnetic field 
spectroscopy, e.g. [116^ 1545]. This will, of course, require very large- aperture installations as the 
relevant energy range is just in the GZK suppression region. 



7 Importance of Modeling Hadronic Interactions 

There are strong indications for shortcomings in the shower simulations, probably due to limi- 
tations of modeling hadronic interactions. 

Detailed studies of the shower development in the atmosphere have been performed with 
the KASCADE multi-detector set-up and interaction models have been improved [214^1208^1546^ 
1547] . A particularly valuable tool to test high-energy interaction models are correlations between 
different shower components j 553| 1554] . Some years ago several models like SIBYLL1.6, DP- 
MJET2.5, or NEXUS 2 failed to describe the measurements of particular correlations. On the 
other hand, for modern models like QGSJETOl, SIBYLL2.1, or DPMJET 2.55, the KASCADE 
measurements are compatible with predictions for various correlations between the electromag- 
netic, muonic, and hadronic components, i.e. the measurements are bracketed by the extreme 
assumptions of primary protons and iron nuclei [553^ I554j . In previous analyses pure proton 
or iron compositions have been assumed as extreme cases. However, at present, more detailed 
analyses are performed [554^ 1555] . They take into account the spectra for elemental groups as 
obtained from investigations of the electromagnetic and muonic components (as discussed above. 



see Fig. 15). These investigations reveal deviations between measurements and simulations for 
the hadronic component of the order of 10% to 20% |554| . 

The situation is similar at very high energy. For example, the mean X^ax of the HiRes- 
MIA data is not consistent with the measured muon densities of the same events [290j . The 
conclusions on mass composition from Haverah Park data are different if the time structure of 
the shower front is used instead of the muon yield that determines the rate of inclined showers. 
The Auger data indicate that the energy scale derived from surface detector simulation seems 
to be of the order of 25% higher than that obtained from fluorescence measurements [357J. 

In the foreseeable future soft multi-particle production will not be calculable within QCD. 
Therefore the modeling of cosmic-ray interactions will continue to depend strongly on the input 
from accelerator experiments that is needed to tune phenomenological particle production mod- 
els. Measurements at both, fixed target and collider experiments can substantially contribute to 
reducing the uncertainties of the models and, hence, determine the composition of cosmic rays. 

High-energy interactions are very important for the overall shower profile but hadronic multi- 



particle production is least known in this energy range. This is illustrated in Fig.|28]in which pre- 
dictions of hadronic interaction models are shown for proton-proton collisions at ^/s = 14 TeV. 
The acceptance ranges for different detector components of the CERN CMS [556j detector for 
the LHC |557| are also shown. 

Every ultra high-energy air shower contains many sub-showers of lower energy. For example, 
the slope of the lateral distribution of particles in a shower is a measure of the mass of the primary 



51 




pseudorapidity r| pseudorapidity r| 

Figure 28: Prediction of the transverse {left) and total {right) energy flow produced 
in proton-proton collisions at LHC as obtained with cosmic-ray interaction models. 
Also shown are the acceptance ranges of the CMS central calorimeter, the CASTOR 
calorimeter and the zero degree calorimeter for neutral particles (ZDC). 

particle. Even for high-energy showers of 10^^ eV, this slope is very sensitive to assumptions on 
hadronic multiparticle production in the energy range of ten to a few hundred GeV jl77| . The 
energy distribution of hadronic interactions, in which at least one meson was produced that in 
turn decayed to a muon that reached sea level, has a maximum in the range between 80 and 
200 GeV. Most of the interactions are induced by pions (70%) and nucleons (20%). 

In addition to accelerator measurements of hadronic multiparticle production [558J, mea- 
surements and understanding of air shower data at lower energy are very important to tune and 
validate the used hadronic interaction models [560J . 

8 Conclusions and Outlook 

The all-particle flux of cosmic rays is reasonably well known up to the highest energies. Recent 
measurements by the HiRes and Auger Collaborations established a GZK-hke suppression of 
the flux at energies exceeding 6 • 10^^ eV. 

In the knee region the mean mass of cosmic rays is found to increase as function of energy. 
The knee is caused by sequential breaks in the spectra of individual elements, starting with 
the light elements. At present, a rigidity dependence of the cut-off energies for the individual 
elements is likely but not yet clear beyond doubt. Above 10^^ eV the situation becomes very 
uncertain. Almost no data are available in the important energy range where the galactic 
cosmic-ray component is expected to end (10^^ — 10^^ eV), at present there are only limited 
experimental efforts in this region. At the highest energies (^10^^ eV), several experiments 
indicate a light to mixed composition with a strong dependence on the model used to describe 
high-energy interactions in the atmosphere. 

Large scale anisotropics have been found at low energies, being compatible with the move- 
ment of the Earth around the Sun (Compton Getting effect). In the knee region, the anisotropics 
disappear, indicating that the rest frame of galactic cosmic rays co-rotates with the Galaxy. At 
ultra-high energies anisotropy measurements provide independent information on the composi- 
tion due to deflections in the galactic magnetic fields. The results of the Pierre Auger Obser- 
vatory indicate a correlation between the arrival direction of cosmic rays and the super galactic 
plane. 
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The energy of the transition from galactic to extragalactic particles is discussed contro- 
versially. The scenarios considered limit the transition energy to the range between 10^^ and 
10^^ eV. Precise composition measurements in this region will be decisive. 

The knowledge about the cosmic-ray composition is presently limited by the uncertainties in 
the hadronic interaction models used to describe the air shower development. The ambiguities 
can not be resolved by cosmic-ray measurements solely. Collaboration with experiments at the 
LHC and fixed target experiments is mandatory to improve the understanding of multiparticle 
production in the (extreme) forward region, thus providing reliable simulation codes. 

Independent and complementary information about the origin of high-energy cosmic rays 
can be obtained by measurements of high-energy gamma rays and neutrinos (multi messenger 
approach). Recent measurements of TeV gamma rays from supernova remnants give strong hints 
for this object class as sources of galactic cosmic rays. Measurements of neutrino and photon 
fluxes are very useful to distinguish between different scenarios for the transition from galactic 
to extragalactic cosmic rays. They are important to establish a GZK feature beyond doubt and 
provide composition information at the highest energies. 

In the near future new cosmic-ray detectors will provide additional information on the transi- 
tion from galactic to extragalactic cosmic rays (KASCADE-Grande [561j . IceTop/IceCube [562j, 
extensions of Auger-South [563! 1564 ], extension of the Telescope Array j565j ), the anisotropy of 
the arrival directions and composition of cosmic rays at the highest energies on the whole sky 
(Auger North [566]) as well as super GZK particles (EUSO [567J, a recovery of the cosmic-ray 
flux beyond the GZK resonance). A promising new detection technique, the measurement of 
radio emission from air showers, is presently developed in the Pierre Auger Observatory and 
the LOFAR project. This technique is expected to allow efficient cosmic-ray measurements in 
future large-aperture experiments. 

In foreseeable future no significant improvement of hadronic interaction models is expected. 
Therefore, the range of direct measurements should be extended to energies approaching the 
knee. Ideal would be an ACCESS-type |569j experiment in outer space with exposure times 
of several years. More experimental efforts are needed in the energy region where the galactic 
component is expected to end. At ultra-high energies the anisotropy studies are limited by 
statistics only. Therefore, large aperture experiments are needed with full-sky coverage such as 
the Pierre Auger Project (with observatories in the southern and northern hemisphere). 
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